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Summary 


In the frequency spectrum of the noise of a fluorescent lamp ballast, “Shum” and “rustle” can 
be distinguished. Hum comprises small even multiples of the mains frequency, at which the 
dimensions of the ballast are small in comparison with the wavelength of the sound in the air. 
The noise at these frequencies is, therefore, mainly produced by the ballast acting as a vibration 
exciter to larger mounting surfaces. Rustle frequencies (about 1000---3000 c/s) are radiated directly 
by the ballast itself. Hum may be eliminated by use of a core construction where the displace- 
ments due to magnetic attraction across the air-gap cancel out the magnetostrictive displacements 
at the mounting points of the reactor. Rustle radiation must be isolated by putting the ballast in 
a closed box. A filling for the empty space in these boxes providing acoustical impact sound 
insulation and good thermal conductivity is described. 


Sommaire 

On peut déceler, dans le spectre des fréquences du bruit produit par le régulateur d’une lampe 
a fluorescence, un «ronflement» et un «bruissement». Le ronflement se compose des petits multi- 
ples pairs de la fréquence du réseau, pour lesquels les dimensions du régulateur sont faibles devant 
la longueur d’onde du son dans lair. Le bruit 4 ces fréquences est donc surtout da au régulateur 
qui fait entrer en vibration les grandes surfaces du bati. Les fréquences de bruissement (environ 
1000---3000 Hz) sont rayonnées directement par le régulateur lui-méme. On peut supprimer le 
ronflement en employant un type de réalisation du noyau tel que les déplacements produits par 
Vattraction magnétique dans l’entrefer soient opposés aux déplacements dus a la magnétostriction 
aux points de fixation de la bobine. On peut empécher le rayonnement du bruit de bruissement en 
plagant le régulateur dans un carter fermé. On indique en quel matériau remplir le restant du 
carter pour assurer un isolement acoustique aux sons d’impact et avoir une bonne conductibilité 
thermique. 


Zusammenfassung 

Im Gerauschspektrum der Drosselspule einer Leuchtstofflampe lassen sich deutlich ,,Brumm*‘- 
und ,,Rausch‘‘-Anteile unterscheiden. Das Brummen umfaft wenige gerade Harmonische der 
Netzfrequenz, fiir welche die Dimensionen der Drossel klein im Verhiltnis zur Wellenlinge des 
Luftschalls sind, so daB der Schall durch Anregung von gréBeren Montagefliichen erzeugt wird. 
Die Rauschfrequenzen (etwa 1000---3000 Hz) werden direkt von der Drosselspule abgestrahlt. 
Das Brummen kann durch Verwendung einer Magnetkernkonstruktion vermieden werden, bei der 
die durch magnetische Krifte iiber den Luftspalt erzeugte Deformation gerade durch die zwischen 
den Montagestellen auftretende magnetostriktive Deformation aufgehoben wird. Die Abstrahlung 
des Rauschens mui dadurch vermieden werden, dafs man die ganze Drosselspule in einen ge- 
schlossenen Kasten setzt. Es wird ein Fiillmaterial fiir den Kasten angegeben, das eine gute 
Schallisolation und zugleich hohe thermische Leitfihigkeit gewahrleistet. 


1. Introduction 


lower [1]. 


makers of fluorescent tube ballasts this 
factory level” is gradually being set lower and 


‘satis- 


some fifteen years ago the problem of the noise 
produced by the associated ballasts was not taken 
very seriously. This new lighting technique was 
used mainly in factories, etc., where the ballast 
noise was usually imperceptible with the general 
sound level obtaining. 

Gradually, however, the use of this type of 
lighting was extended to offices, shops and homes. 
These quieter surroundings made a reduction of 
the noise below a “satisfactory level” imperative. 
As a result of competition between the various 


In this paper a survey will be given of some 
techniques of noise abatement in fluorescent 
lamp ballasts, resulting from the research work 
done on this subject at the N. Y. Philips’ Gloei- 
lampenfabrieken. 


2. Character of the sound 


The sources of sound in a fluorescent lighting 
installation are the ballasts. These contain iron- 
core inductors, usually with air-gaps in the 
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magnetic circuits. The magnetic field in an iron 
core produces mechanical deformations of the 
latter, both by the action of attraction forces 
across the air-gap and by magnetostrictive 
effects. The periodic character of the magnetic 
field makes these deformations result in vibra-' 
tions and noise. 

The noise as produced by a fluorescent light- 
ing installation is distributed over a large part 
of the audible frequency range. We shall sub- 
divide it into “hum”, “rustle’”’ and (occasionally) 
“rattle’’. 

Rattle is an indirect consequence of vibrations 
in fixtures, etc., of which many parts are held 
together by the force of gravity only. As soon as 
the peak value of the acceleration amplitude of 
vibrations in the construction exceeds the gravi- 
tational acceleration, those parts may start ratt- 
ling against each other. Such vibrations are often 
caused by impacts on the ceiling on which the 
fixtures are suspended. Older types of ballasts 
occasionally produced, during operation, rather 
violent vibrations with the same effect, but with 
modern ballasts there is little chance of vibrations 
being excited which come anywhere near the crit- 
ical level for producing rattle. We shall there- 
fore consider rattle to be outside the scope of 
this paper. 

In the remaining sound two important fre- 
quency regions may be discerned, leading to the 
division into “hum” and “‘rustle’’. 

Hum concerns the low-frequency region up to 
a few hundred c/s within which discrete fre- 
quencies occur which are even multiples of the 
mains frequency. In comparison with the wave- 
length of these tones in the air a normal ballast 
is so small that, if suspended elastically, it will 
radiate hum frequencies as a point source in 
accordance with its total volume displacement. 
This radiation is negligible in comparison with 
the sound production obtained if the ballast acts 
as a vibration exciter to larger surfaces (walls, 
ceilings, large parts of fixtures, etc.) which, as a 
consequence of their dimensions not being small 
in comparison to the wavelength, are capable of 
radiating the hum noises. 

Rustle indicates the sound within a frequency 
range of about 1 to 3 ke/s where, in the sound 
analyses, rather surprisingly high intensity 
levels may be found. Unlike the low frequency 
vibrations these high rustle frequencies are 
radiated to an appreciable extent by the ballast 
itself, as the source is no longer small compared 
with the wavelength. Moreover, the frequency 
range of the rustle noise is such that natural 
mechanical vibrations of the coil construction— 
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although heavily damped—may add to the irreg- 
ular character of the phenomenon. 

Measurements of the sound radiated by an 
elastically suspended ballast, such as have been 
made by Mason and Nutr [2], will reveal only 
its rustle properties. Only if the dimensions of the 
reactor are large enough can this method be used 
to study hum. Large transformers as used in 
distribution networks of electrical energy have 
been studied in this way by Gorpy [3]. GouLp 
[4] tried to overcome the difficulty by studying 
the noise of fluorescent-lamp ballasts when 
mounted on a rather large fixture, thus assuring 
at least some hum radiation. We have investi- 
gated the problem of.hum by studying the low- 
frequency vibrations of the iron cores of the 
ballasts since the only effective way to reduce 
hum and avoid rattle is the reduction of vibra- 
tion transmission from the ballast to the base on 
which it is mounted. In addition, rustle frequen- 
cies, in so far as they are radiated by the ballasts 
themselves, must be isolated. 

It will be understood that the principal aim of 
this research was to find measures for noise reduc- 
tion which can be taken in the manufacture of 
the ballasts. These measures therefore can only 
be those of such a nature as will give satisfact- 
ory reproducible results without the need of in- 
dividual acoustic adjustments of the ballasts. (This 
would become too costly for mass production.) 


3. Theory of hum reduction 


To study the action of a ballast as a vibration 
exciter, especially at low frequencies, we shall 
analyse the causes of mechanical vibrations in a 
magnetic core. 

If an iron core, generally comprising an air-gap, 
is magnetised two deformations oceur simul- 
taneously, one elastic and the other magneto- 
strictive. 

The elastic deformation originates in the force 
with which the limbs of the iron circuit attract 
each other across the air-gap. This force F in 
newtons is given by 
107 B?O 


= 
82 


(N), (1) 
where B = the magnetic induction in Vs/m* and 
O = the cross section of the magnetic circuit in 
m?, For an induction B varying sinusoidally with 
time: B= Bnaxsinwt, and the force F' becomes 
proportional to (1+ cos2mt). This means that 
the magnetic attraction varies sinusoidally bet- 
ween zero and twice the average value, with a 
frequency twice that of the current in the coil. 
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Hum-frequencies being low in comparison with 
any natural frequencies of the construction, the 
momentary deflection pattern and the moment- 
ary force will be related in accordance with the 
static case. 

The magnetostrictive properties of the magnet- 
ic material cause relative elongations and com- 
pressions of the core in the direction of the 
magnetic flux. The technical aspect of the phen- 
omenon has been studied, amongst others, by 
FAHNOE [5], LEONARD [6] and Wixu1AMs [7]. The 
characteristic relation between the specific mag- 
netostrictive elongation ¢«=AL/L and the induc- 
tion B for ordinary transformer iron is given in 
_ Fig. 1. At values of B up to about 80 % of satura- 
tion, € is positive (magnetostrictive elongation) 
and may be taken as approximately proportional 
to the energy density of the magnetic flux: 


107 B2 
oo ps (2) 
where « = a proportionality factor. 


For ordinary transformer iron, depending upon 
the Si-content and the treatment, « will be of the 
order of magnitude (1 to 3)-10°1 m?/N. 

Eq. (2) holds within the range wherein the 
wrongly directed magnetisations within individ- 
ual Weiss domains are reversed. At high induct- 
ances the magnetisation within Weiss domains 
changes into less preferred crystallographic direc- 
tions. This produces a magnetostrictive shorten- 
ing which soon supersedes the elongation already 
existing till, at the saturation value, ¢« amounts 
to about —7-10-®. In normal operation the 
induction in a core will seldom exceed the range 


wherein Eq. (2) holds. 


+310° 


| Bsat 
| 0 


10° 


Fig. 1. Characteristic relationship between magnetostric- 
tive elongation ¢ and induction B. 
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Eq. (2) indicates that, like the elastic deforma- 
tion, the magnetostrictive elongations will vary 
sinusoidally at twice the frequency of the (sinus- 
oidal) coil current. As the magnetic design of 
inductors generally aims at a uniform flux 
density, the magnetostrictive deformation pat- 
tern will to a first approximation be that of a 
proportional increase of all dimensions in a plane 
parallel to the core laminations, without bending 
effects. 

Actually some bending due to magnetostrictive effects 

is unavoidable, since at the corners of the core the flux 

density at the inside of the bend will be greater that at 

the outside. (Cf. Stern, [8]). As a result, for instance, a 

normally right angle between the limbs of a magnetic 

core will be a little more than 90° at the moments of 
maximum induction. This leads to some bending of the 
straight parts of the core. 

We shall now consider a core construction, a 
sectional view of which is given in Fig. 2. The 
air-gap is situated symmetrically at the centre 
of gravity of the whole. The core is considered to 
be elastically suspended, so that it may move 
freely without exerting any external forces. The 
centre of gravity then remains stationary. 


3, 
I 
ee NY 


: 
rx 
i 


| MM OS 


Fig. 2. Sectional view of symmetric core construction. 


If the core is magnetised, the magnetostrictive 
displacement of points along the line P—P’ will 
increase in an outward direction at both sides of 
the air-gap, where it is zero, due to the constant 
specific elongation of the centre limbs (the flux 
density along which is assumed to be constant). 

The elastic displacement of points along P—P’, 
on the contrary, is largest in the planes bounding 
the air-gap, and decreases in outward directions 
due to the elastic elongation of the centre limbs. 
As a consequence, somewhere along the line 
P—P’, two points exist—symmetrically situated 
with respect to the air-gap—where the total 
displacement due to both effects is zero. These 
points may lie outside the geometrical boundaries 
of the core. By a proper adjustment of all para- 
meters these points can be made to coincide with 
the points P and P’, where the core can then be 
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clamped to its mounting device without any 
vibration being transmitted from the core to the 
base. 

The following symbols will be used: 


6 (m) = length of air-gap in the direction of: 
the flux, 

O(m?) == cross-section of the iron core at the 
air-gap, 

O'(m?) = cross-section of the non-magnetic 


material in the gap (O’ < 0), 

L (m) = distance between the fastening 
points P and P’, 

E.(N/m?) = Young’s modulus of the magnetic 
iron, 

E, (N/m?) = Young’s modulus of the non-mag- 
netic material in the gap, 

R.(N/m) = rigidity of the iron core at the place 
of the gap, against elastic deforma- 
tion in the direction of the mag- 
netic flux. 

The total rigidity R1o: of the iron core together 
with the gap filling material, against elastic de- 
formation in the direction of the magnetic flux 
becomes: ; 

Root ae R. os a 
0) 

Under the action of the attractive force bet- 
ween the limbs bounding the air-gap, as given 
by Eq.(1), the elastic deformation of 6, Ad, 


amounts to: 


Aé = 


(N/m). (3) 


UDO) I 
Sa Ree ey 


The negative sign indicates a contraction of 
6. This contraction corresponds to a force F, = 
—R. Ao in the iron core and a force F, in the gap 
material (F = F,+ F,). The force F. produces 
an elastic elongation (AL)’ of the centre limbs 
bounding the air-gap. Supposing LS6 this 


elongation is given by: 


aC: L 
(AL)'= —R. Ad E.0 (m). (5) 
The total elastic elongation of L, (AL)ectast. 
= Ad+ (AL)’, hence with Eq. (4) and (5): 
LOTR a 1 R.L 
Ba \R.+E,078)\1~ £0) ™- 
The magnetostrictive elongation of the distance 
P—P’, (AL) magn. according to Eq. (2), amounts to: 
107 B? 
(AL) magn. ——W F (oa ay a (m). (7) 
From (6) and (7) the total elongation AL of the 
distance P—P’ becomes: 
(8) 


107B2 0 RL 
War a x [bs feet ro) | 


(AL) ctast. or 
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There are three important particular cases to be 
considered : 

Case 1: 6=0 (no air-gap; purely magneto- 
strictive deformation). From (8), with 6 =0 we 
see: 


(AL) 5—0 = -La (m). (9) 

This formula can be used for determining the 
magnetostrictive constant « of the type of core 
material to be used in the production of the 
ballasts. A special core is made of the intended 
shape, but without an air-gap. The deformation 
of the easily accessible points P and P’ is meas- 
ured as a function ofthe induction B, from which 
x may be calculated. 

Case 2: E, = 0 (no filling material at all in the 
air-gap). In this case: 

7 B2 

(AL)z,—0 = ae (toe 5 £) (m). (10) 


With the numerical value of « known from 
case 1, and of FE. as derived from tables of elastic 
constants, measurements of AL as a function of 
B on a special core without air-gap filling will 
produce a numerical value of R.. The core for 
this purpose must be identical in shape with the 
production type, with the filling material in the 
gap omitted. 

Case 3: AL =O (the ideal production type). 
In this case Eq. (8) leads to: 


Ej00cyEeO Sah 
re) aLek, 


R. (N/m). (11) 
With the numerical values of « and R, as found 
in the cases 1 and 2 known, the value of E, to 
fulfil the condition of no vibration at P and P’ 
can now be calculated. 

The method of assuming the value of L and 
calculating the corresponding value of E, has 
been chosen because in practice the clamping 
points of a given core construction are more or 
less fixed, whereas the value of E, may be varied 
to a much wider extent. 


Numerical example 

For a given type of core stampings, to be made 
into ballasts for 40 W fluorescent tubes on 225 V, 
50 c/s mains the dimensions (cf. Fig. 2) were: 
a= 108 mm, b=58mm, c=17mm,6=—1.3mm, 
0 = 17-23 = 390 mm’?, L= 97 mm, E.=2-104 
N/m?. 

Measurements on a core without air-gap gave 
(AL) magn. = 0.108 um at B=1.25 Vs/m?, hence 
with Eq. (2): «= 1.8-10-! m?/N. 

Similar measurements on a core without filling 


material in the gap gave (AL)z,-o = —0.096 ym 
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at B=1.25 Vs/m?, from which, using Eq. (10), 
we may calculate R.=4.8-108 N/m. 

On inserting these values into Eq. (11) we find 
that the ideal material to fill the entire gap 
(O’= O) must have a value of Young’s modulus 
E,=1.4-10° N/m?. A choice of the proper filling 
material for the gap of the intended core can 
thus be made. 


4, Measuring apparatus 


The apparatus used for measuring AL is shown 
schematically in Fig. 3. A pair of electrodynami- 
cal vibration pick-ups of the type “Philips” 
PR 9261 [9] are used in u-metal boxes for 
screening against the leakage fields of the ballasts 
under test. 

The moving parts of these vibration pick-ups 
possess a mass of about 10 g each and are pressed 
against the iron core with a constant force of 
about 8N, thus the maximum measurable 
acceleration becomes 800 m/s’. By using two 
pick-ups the constant pressing forces oppose each 
other. The inductor to be measured is suspended 
on rubber springs and can thus move freely, the 
only additional external forces being those exerted 
by the probes of the vibration pick-ups. The 
resulting constant compression of the core 
construction, although sometimes of the same 
order of magnitude as the measured deformation, 
is of little importance as the measured sinusoidal 
variations of the core dimensions are caused only 
by the alternating magnetisation. 


PR 9250 
with 100 c/s 
calibration voltage 


Selective 


amplifier 


Fig. 3. Schematic arrangement of measuring apparatus. 
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In many cases the pick-ups can touch the core 
at symmetrically situated points (cf. points P 
and P’ in Fig. 2), which, during operation of the 
freely suspended coil, perform vibrations of iden- 
tical amplitudes and phase, in opposite directions. 
In these cases the voltages of both pick-ups, indi- 
cating the velocity amplitudes at the measured 
points, are identical in phase and amplitude, and 
may be combined electrically, as shown in Fig. 3, 
in a network where V= 0.5 (V,+ V,). This is a 
special feature to increase the workability of the 
set-up, since undesired vibrations in the meas- 
uring locality cause the suspended coil to act as 
a pendulum, thus producing identical voltages in 
opposite phase in the two pick-ups. These are 
cancelled out in the averaging network, so that 
the output voltage of this network indicates the 
velocity amplitude of the “breathing” of the 
core, with the background disturbances very 
much reduced. In order to obtain the actual core 
deformations, this voltage is integrated electri- 
cally in the calibration unit PR 9250 [9], [10], 
and after amplification in a selective amplifier, is 
fed through one channel of the electronic switch 
to the cathode-ray oscilloscope (Fig. 3). The 
selective amplifier consists of two successive 
stages with negative feedback through R-C 
filters of the bridged-T type, and passes a narrow 
band at 100 c/s. The other channel of the elec- 
tronic switch is used for simultaneously indi- 
cating the current in the measured coil. Thus one 
can see whether an expansion or a contraction 
occurs at the moments of maximum induction. 

In the normal use of a PR 9261 vibration pick- 
up in conjunction with the calibration unit 
PR 9250, the amplitude of the measured vibra- 
tion is found by comparing the amplitude on the 
oscilloscope with that due to a 50 c/s comparison 
voltage supplied by, and read on, the calibration 
unit. Because of the selective amplifier in the 
circuit, our comparison voltage has to be a 100c¢/s 
voltage, or rather at all times exactly twice the 
mains frequency because of the steep slopes of 
the pass band in the selective amplifier. The 
calibration unit was therefore supplied with a 
frequency doubler, together with an attenuator 
so that the actual calibration voltage at 100 c/s 
was exactly one hundredth of the value indi- 
cated on the scale of the instrument. 

With the apparatus described a sensitivity of 
3-10-°m of the vibration amplitude measure- 
ments was obtained for symmetrical cores. In 
non-symmetrical cases, where the voltages of 
each pick-up must be measured separately, the 
background disturbances determine the sensitiv- 
ity. In this case one has to work during very 
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quiet periods, e. g. at night, to obtain anything 
like the sensitivity value mentioned above for 
symmetrical cases. 


5. Less symmetrical core constructions 


The example of treatment for the minimizing 
of hum vibrations given earlier was simplified in 
that the shape of the core (Fig. 2) was symmetri- 
cal with respect to the line P—P’. Often, for 
reasons of economy in mass production or other- 
wise, such a symmetrical design is not feasible. 
An example of a less symmetrical construction 
encountered in practice is given in Fig. 4a. 

The only line of symmetry in this case is Q—Q’ 
through the central air-gap. If P and P’ again be 
the clamped points we now have to consider 
separately the components of vibration of P and 
P’ along the line P—P’ and perpendicular to this 
line. In view of the symmetry with respect to 
Q-Q’ the treatment of the vibration components 
along P—P’ is analogous to that described before. 
Eq. (8) may be used and thus the ideal value of 
Young’s modulus for the filling material in the 
air-gap may be determined. As would be expected, 
the value of R. proves to be an order of magni- 
tude inferior to that in analogous constructions 
of the symmetrical type pictured in Fig. 2. 

For example an inductance to be used in series 
with a 40 W fluorescent tube on a 225 V, 50 c/s 
mains possessed the following dimensions (cf. 
Fig. 4a): a = 108 mm, b = 45 mm, c = 17 mm, 
6=13mm, L=97mm, 0= 17-23 = 390 mm?2, 
Ke 2 10 IN ime: 
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Fig. 4a. (above) Deformation pattern of non-symmetrical 
core with symmetrical gap-filling. 
Fig. 4b. (below) Ditto with non-symmetrical gap-filling. 
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Measurements on a gapless core gave ~ = 
27-10-13 m?/N; measurements on a core without 
gap filling showed a component of AL along P—P’ 
at B=1.25 Vs/m? of —6.6um (contraction), 
from which may be calculated (according to Eq. 
(10)): R.=3.4-10’ N/m. 

The correct material to fill the whole cross- 
section of the air-gap (O’= O) must then possess 
a value of Young’s modulus E, = 4.6-10° N/m?. 

An inductance equipped with this air-gap 
filling will still show vibrations of the points P 
and P‘ in a direction perpendicular to P—P’. This 
is to be expected, since the back B—B’ of the core 
will expand due to magnetostriction, whereas the 
attractive foree across the air-gap, apart from a 
small contraction, will exert mainly a bending 
torque along B-B’. Bearing in mind that the 
distance between P and P’ will not alter during 
magnetisation, the centre lines of the core limbs 
can be expected to be deflected as shown (exag- 
gerated) by the dotted line in Fig. 4a. The shape 
of this dotted line was derived from measurements 
on a ballast the coil of which was divided into 
two parts, one on each of the limbs R-S and 
R’-S’, so that the core was accessible for the 
vibration pick-up probe near the air-gap, around 
the sides and at the back. The deformation com- 
ponents were measured at a number of points in 
two perpendicular directions, with B= 1.25 
Vs/m? and the resultants drawn in vectorially. 
The scale of the deformation is also given in 
Fig. 4. 

The shape of the deflection pattern indicates 
that the distance between P—P’ and the centre 
of gravity of the construction (situated some- 
where along Q—Q’) will vary, or in other words, 
for a freely suspended ballast of this type, the 
clamping points P and P’ vibrate perpendicularly 
to the line P—P’, since in this case the centre of 
gravity of the whole construction remains at rest. 
Unfortunately, due to the vibrational movements 
of the core, this centre of gravity is not a fixed 
point. Moreover, the total mass of such an indue- 
tor is about equally distributed between the iron 
core and the coil, hence any mutual displacement 
of these two will have an appreciable effect upon 
the exact momentary position of the centre of 
gravity of the whole. 

The coil in itself is a rather stiff tube which is 
clamped around the limbs R-S and R’-S’ with 
certain mechanical tolerances. The vibration 
amplitudes involved are small with regard to 
these tolerances, and it is therefore purely a 
matter of chance whether the clamping between 
core and coil will be most effective at R and R’, 
or S and S’, or somewhere in between. As a result 
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the mutual displacements of coil and iron core, 
and therefore also the amplitude of the remaining 
vibrations of P and P’ perpendicular to the line 
P-P’ will be unpredictable and, moreover, will 
vary considerably amongst ballasts in a given 
production batch. 


The shape of the dotted line in Fig. 4a indicates 
that a possible remedy against mutual displace- 
ments of core and coil would be a bending torque, 
_ proportional to the acting attractive forces, and 
applied in such a way that the limbs R-S and 
R’-S’ would stay aligned during deformation. 
This may be obtained, for instance, by filling the 
upper half of the air-gap (cf. Fig. 4b) with a 
material of which E, is lower than the calculated 
value, and the lower half of the gap with a 
material with a high value of E, so that the 
overall mechanical rigidity of the gap filling is 
_ the same as if the whole gap were filled with the 
correct material to fulfill the condition AL = 0. 
Now, however, the seat of the resultant of the 
reactive force opposing compression of the gap 
filling will be shifted towards the side where the 
hardest filling material is located, and this force 
together with the magnetic attraction, acting in 
the centre line of the magnetic circuit, may pro- 
vide the desired bending torque. 


By varying the method of gap filling, always 
fulfilling the condition AL = 0 for vibrations 
along P—P’, it is thus possible to bring the reac- 
tive force of the gap filling out of centre to a 
variable extent, and it is thus even possible to 
make cores the limbs R-S and R’-S’ of which 
bend outwards during deformation. For a given 
type of core construction the right combination 
of materials in the gap, which both makes AL = 0 
and keeps the coil limbs aligned, may easily be 
found by measurements of the deformation 


pattern (Fig. 4b). 


For a ballast where these conditions are ful- 
filled the coil will not move perpendicularly to 
the line P—P’, no matter where it is clamped. 
Nevertheless P and P’ will still show some per- 
pendicular vibration with respect to the centre 
of gravity of the whole, due to the unilateral 
bending of the back B—B’ of the magnetic circuit 
(which still exists, notwithstanding the bending 
torque). A little over-compensation of the coil 
movements, to be obtained by making the limbs 
R-S and R’-S’ bend outwards, would appear to 
be attractive. This, however, again introduces 
variations due to the small differences in the 
location of the most effective clamping point for 
the individual coils in a production run. There- 
fore a good alignment of the limbs, with constant 
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results as to the hum properties of the ballast, 
proved to be acceptable. 


6. Transfermers 


A fluorescent-lamp ballast for use on low- 
voltage mains (110 or 125 V) consists of a step-up 
transformer with leakage path for the magnetic 
flux, thus combining, in effect, a transformer and 
a choke. Occasionally small transformers are 
used separately. For instance fluorescent lamps 
of the Philips’ type “TL” S, which operate on 
225 V mains in series with an incandescent lamp 
as ballast, require a step-up transformer on 110 
or 127 V mains owing to the higher starting 
voltage needed. In incandescent spotlight in- 
stallations for shop window lighting step-down 
transformers are sometimes used of a size com- 
parable with fluorescent lamp ballasts. 

Transformers do not possess an air-gap and 
thus the technique of compensating the elastic and 
magnetostrictive deformations can not be applied. 
In operation a transformer core shows only mag- 
netostrictive deformation. To avoid the trans- 
mission of vibrations to the mounting surface we 
can either use a transformer iron with little or 
no magnetostrictive effect, or brace the trans- 
former exactly at its centre of gravity. Both these 
measures are, however, impracticable. We must 
therefore consider the fact that all points that 
may be used as clamping points vibrate during 
operation. 

Mounting devices that tend towards a freely 
suspended state of the transformer are out of the 
question. If they are very supple, so as not to 
transmit the vibrations, they are unable to hold 
the transformer firmly in place, which is a basic 
practical condition in the layout of an installa- 
tion. 

However, if two points exist on the core which 
can be used as clamping points and which, in a 
freely suspended state of the transformer, vibrate 
in the same line with identical amplitudes and 
phases but in opposite directions, it is possible 
to design a mechanical device providing a firm 
fastening of the whole construction whilst pre- 
venting the transmission of magnetostrictive 
vibrations towards the mounting surface. 

Two types of fastening members have been 
developed for this purpose. They are the U-shaped 
bracket with converging legs and the rectangular 
U-shaped bracket with hinged legs. 

The convergent-U-shaped bracket is pictured 
schematically in Fig.5. Only the elastically 
neutral lines of the legs and the body are drawn. 
The transformer is supposed to be of the type 
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Fig. 5. Sectional view showing deformation of elastically- 
neutral lines of convergent-U-shaped bracket. 


illustrated in Fig. 2 but without an air-gap. The 
clamping points P and P’ are braced by the ends 
of the legs of two convergent U brackets, one on 
either side of the iron core. The centres F of the 
bodies of both brackets are fastened onto the 
mounting surface. 

At the moments of maximum magnetic induc- 
tion B the points P and P’ are displaced out- 
wards over a distance 0.5 AL each. If the clamp- 
ing is firm enough the ends of the legs of the 
brackets will move with the points P and P’ and 
parallel to themselves. The dashed line in Fig. 5 
will then indicate schematically, and in an exag- 
gerated way, the shape of the elastically neutral 
line of the bracket in the deformed position. 

Calculation of this elastic deformation leads to 
the following expression: (12) 


AH =: (L/H)? cosy + (L/H) siny — siny - tany — (J,/I,) tany 
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Fig. 6. Graph of values satisfying Eq. (13). 


Eq. (13) gives the ratio between the moments 
of inertia of the body and the limbs of the bracket 
required for absence of vibration transmission 
from the transformer to the mounting surface. 
In Fig. 6 this ratio is plotted graphically against 
L/H with y as parameter. The equation holds if 
the mechanical resonance frequencies of the 
construction are high with respect to the vibra- 
tion frequency concerned (mostly 


100 c/s). This can easily be arrang- 


ae 

where 

L = the distance between P and P’ (m) (cf.Fig.5) 

H = the distance between the line P—P’ and the 
elastically neutral line at the centre F of 
the body of the bracket (m), 

AL and AH the increases of L and H, respec- 
tively, during deformation (m), 

y = half the angle between the converging 
limbs, 

I, and I, the moments of inertia to bending in 
the plane of drawing of Fig. 5, resp. of the 
body and the limbs of the bracket (m‘). 

AH may thus be positive or negative. Sup- 
posing that the mass of the bracket is negligible 
compared with that of the transformer, the 
centre of gravity of the latter, which is halfway 
between P and P’, will also apply to the whole 
construction. In this case no vibration will be 
transmitted if the distance between F and the 
common centre of gravity is invariable, in other 

words if AH = 0. By inserting this into Eq. (12) 


we obtain: 


4 (L/H) cosy + 4-siny + I,/I, 


ed, as only the ratio between the 

moments of inertia I, and I, en- 
ters into Eq. (13), and not the absolute values. If 
fairly strong members are used for fastening, whilst 
still satisfying Eq. (13), it is possible to obtain a 
construction which is rigid’ enough for practical 
use, and of which the lowest natural mechanical 
resonance frequency is high enough for the 
suppositions underlying Eq. (13) to be valid. A 
bracket of this type is pictured in Fig. 7a. 

The rectangular bracket with hinged legs is 
pictured in Fig. 7b. In comparison with Fig. 7a 
it has a very simple appearance, yet may proye 
satisfactory as to lack of vibration transmission. 
Its design is based on the characteristic deforma- 
tion pattern of a symmetrical transformer core. 
Both the body and the legs of the bracket are 
stiffened by flanges, but the corners are inten- 
tionally kept supple, so that they can be regarded 
as hinges in comparison to the stiffened members. 

The deformation, as measured at B= 1.25 
Vs/m*, of the transformer core for which the 
bracket has been developed, is pictured in Fig. 8. 
The bending of the sides, a reason for which 
has been given earlier [8], is also shown in the 
figure. 
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Fig. 7a. (above) Transformer equipped with convergent-U- 
shaped bracket. 

Fig. 7b. (below) Transformer equipped with bracket with 
hinged legs. 
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Fig. 8. Sectional view 
showing defor- 
mation pattern 
of transformer 
core. 
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The centre parts of the ends, between X and Y, 
and X’ and Y’ respectively, show little bending 
but move practically parallel to themselves. If 
the legs of the bracket bracing the core at these 
places also move parallel to themselves, the body 
of the bracket will be stretched without bending, 
so that its centre F remains stationary. This 
centre is again fastened to the mounting surface. 
The forces required for stretching the body 
convey bending torques to the clamped areas of 
the legs, which will prevent them from moving 
parallel to themselves unless the fastening is so 


E.W. VAN HEUVEN: NOISE OF FLUORESCENT LIGHTING INSTALLATIONS 


109 


tight that the clamped areas of the legs and the 
core are mechanically a whole. This can be achiev- 
ed with quite a moderate degree of friction 
between the parts, as can be shown both by cal- 
culation and by experiment. 


7. Rustle abatement; canned ballasts 


We have dealt up to now with measures to be 
taken against the transmission of vibrations from 
the ballast to the mounting surface. This is suffi- 
cient to eliminate the hum noise of the installa- 
tion, but not the rustle, i.e. the less important 
part of the noise, comprising frequencies in the 
range | to 3 ke/s, which are radiated directly as 
sound by the vibrating core itself. 

Rustle originates in the electrical series con- 
nection of an iron cored choke with a gaseous 
discharge. This leads to two phenomena, both 
resulting in core vibrations at rustle frequencies. 

At the moments of re-ignition of the discharge, 
once each half-cycle of the mains voltage, a pulse 
excitation occurs which causes damped mechani- 
cal oscillations of the core at its natural mechani- 
cal resonance frequencies. The lower frequencies 
of these fall approximately in the range of rustle 
frequencies referred to. If, in the measuring 
apparatus, the selective amplifier is replaced by 
a non-selective pre-amplifier so that the rustle 
frequencies are passed, damped oscillations will 
be perceptible on the oscilloscope during about a 
tenth of a half cycle of the mains frequency. The 
higher resonance frequencies are damped more 
heavily, in addition to being excited to a lesser 
extent owing to the finite width of the pulse. 

The other phenomenon can be separated by 
short-circuiting the starter of a fluorescent lamp, 
thus leaving a circuit of two fluorescent lamp 
filaments in series with the ballast and fed from 
a variable, sinusoidal, mains voltage. At low 
values of the latter a pure sinusoidal current 
flows. If the voltage is increased until the peak 
value of the voltage drop across each filament 
exceeds a critical value of about 17 V, discharges 
are set up in the neighbourhood of the filaments 
and at the same time the ballast starts to rustle. 
A study of the current reveals that the discharges 
are random pulses, which are quenched by the 
shunting effect of the filament bringing down the 
voltage drop. 

In normal operation of the fluorescent tube 
this phenomenon may also exist, possibly in 
combination with plasma oscillations. The ran- 
dom pulses of the tube current produce mechani- 
cal core vibrations in the same way as do the 
re-ignition pulses of the main discharge. 
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It is interesting to note that both rustle 
phenomena have a recurrent frequency of 100 
times per second, one being associated with the 
re-ignition of the discharge, the other with the 
peak tube current. An elastically suspended 
ballast, which can practically only radiate its 
rustle as sound, may yet be judged to produce 
“hum”. This sensation of the recurrent frequency 
of the high pitched rustle noise may be another 
evidence of SCHOUTEN’S hypothesis concerning 
the “residue” in sounds with many spectral 
components [11], [12], [13]. 

The strong interaction between tube and 
ballast in the origin of rustle noise is also apparent 
when the sound emitted by an elastically sus- 
pended ballast is measured [2]. When selecting 
ballasts for their rustle properties one may 
compare the sound levels of various ballasts 
feeding successively the same fluorescent lamp, 
but an equally wide variety of sound levels are 
found if different fluorescent lamps are connected 
successively to the same ballast. 

What can the ballast manufacturer do to mini- 
mize rustle noise ? The special clamping methods 
and devices employed to overcome hum excita- 
tion of the mounting surface do not apply in 
general, since the phase relationships of the 
vibrations in the core are substantially different 
from those at frequencies low with respect to the 
natural frequencies of the construction. 

This is, however, of little concern as it is very 
difficult to excite vibrations in building construc- 
tions at 100 c/s or higher. Exxrne [14] found 
attenuations of some 30 dB at the contact point 
even when precautions were taken to assure a 
good contact. 

The main concern is, therefore, the isolation of 
the sound emitted directly. A loose cover over the 
ballast may have some effect but generally leaves 
too many leaks to be a real insulator in the 
acoustical sense of the word. The best thing to do 
is to put the ballast in a closed box (canned 
ballasts). The canning of ballasts presents a 
number of problems, mainly of a thermal nature. 
There must be a good thermal conductivity from 
the ballast to the box for dissipating the deve- 
loped heat without the temperature of the coil 
becoming too high. Metallic contacts between the 
ballast and the box must be avoided as the thin 
metal wall of the latter can be much more easily 
excited at rustle frequencies than building 
structures, in which case its acoustical effect 
would be spoilt. A well-known technique is to 
have the ballast float in a bituminous filling in the 
box. Since bitumen is a poor conductor of heat, 
a mixture of bitumen and sand is often used. 
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The difficulty is that the bitumen must be 
pourable for filling the can, at a temperature not 
detrimental to the electrical insulation. At the 
normal operation temperature of the ballast it 
will therefore tend to become rather soft and a 
really floating ballast will sooner or later come 
into contact with the wall of the can. Additional 
materials are therefore required to keep the 
ballast in its place, but these will increase the 
acoustic contact between ballast and box. 

A good solution of the acoustical and the 
thermal problem has been found by filling the 
empty space in the can with corrugated alumi- 
nium foil. By choice of a suitable “filling factor” 
a volume heat conductivity can be procured that 
equals or surpasses that of bituminous-sand mix- 
tures. Vibrations at rustle frequencies are ade- 
quately attenuated by the succession of many 
contact points between the foils. 

At hum-frequencies the core shows places of 
large and of small amplitudes which all come 
into contact with the aluminium foil initially. If, 
thus, large vibrational forces were to be trans- 
mitted through the foil to the box and its 
mounting surface, it is very likely that Hooke’s 
law would be surpassed locally at the contact 
points of the rather ductile aluminium foil, thus 
providing room for the core vibrations. The foil 
thus tends to hold the ballast, ultimately, at its 
nodal points only. 

In this connection, it has occasionally been 
found that altering the position of a ballast of 
this type with respect to gravity produces a 
slight increase of its hum level, which soon dies 
down to the original value. 


Conclusion 


The noise that will be produced by a given 
ballast is largely dependent upon the cireum- 
stances under which the ballast operates. Even 
the acoustically-best ballasts may produce audible 
hum if properly placed on a tuned organ pipe or 
kettledrum. Equally, these same ballasts may be 
heard to rustle if connected to a fluorescent lamp 
selected for its high rustle output. 

If one wishes to establish a standardized 
quality indication for noise of fluorescent tube 
ballasts in the same way as tests on the electrical 
behaviour have been standardized [15]---[18], 
one should therefore create standard conditions 
under which the ballast is operated during the 
noise test. If these standard conditions comprise 
an elastic suspension within a box of given volume 
and acoustic absorption, as for instance in the 
apparatus used by Mason and Nourtr [2], a 
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sensible indication of rustle may be obtained but 
the problem remains of defining a fluorescent 
tube producing a standard rustle. Dummy 
resistances, as prescribed in the electrical tests 
of ballasts, are of no use in this connection. 

In addition, the vibration transmission at hum 
frequencies should be measured with the ballast 
mounted on and loaded by a specified mechanical 
impedance. Finally, the relative importance of 
these two data must be decided upon, so as to 
obtain a single—weighted— combined figure. 


The author wishes to acknowledge the work of 
many colleagues at N. V. Philips’ Gloeilampen- 
fabrieken, from the results of which this paper 


has been compiled. 
(Received 27th Sept., 1954.) 
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International Acoustics Congress in U.S.A. — 1956 


The Second Congress of the International Commission 
on Acoustics will be held June 17--:24, 1956 in Cambridge, 
Massachusetts, U.S.A., in conjunction with the Fifty-first 
Meeting of the Acoustical Society of America. This is one 
of a series of triennial congresses sponsored by the Com- 
mission [1]. 

The First ICA Congress [2], held in the Netherlands in 
June of 1953, was attended by some 300 scientists and 
engineers from 20 different countries. In addition to tech- 
nical papers and excursions, Congress participants and their 
families enjoyed sightseeing trips, receptions, and other 
functions which formed a part of the Congress program. 

The Commission extended an invitation to the Acousti- 
eal Society of America to organize the 1956 Congress in the 
United States. The Council of the Society, in turn, pro- 
posed that the Congress be held in conjunction with a 
regular meeting of the Society. In this way, all members 
of the Society would have an opportunity to participate in 
the International Congress, and all of the Congress members 
from other countries would automatically join in a meeting 
of the Society. The Council selected Cambridge, Massa- 
chusetts, as the location for the meetings and designated 
Professor R.H.Borr as Chairman of the Organizing 
Committee. 


The central themes of the technical program will be 
Physical Acoustics and Bioacoustics, with special emphasis 
on the inter-relations between these fields and on practical 
applications. Technical sessions will be held at the Massa- 
chusetts Institute of Technology and Harvard University. 
The program will include a number of visits to points of 
special acoustical interest, as well as entertainments and 
opportunities for sightseeing. 

Detailed plans for the Congress will be announced in 
special bulletins. If you would like to attend the Congress 
and wish to receive further information please write to: 


Mr. Joun A. Kesser, Secretary 
Second ICA Congress 

Acoustics Laboratory 

Massachusetts Institute of Technology 
Cambridge 39, Massachusetts, U.S.A. 


{1] The organization and functions of ICA are reported in 
Acustica 1 [1951], 143 and J. acoust. Soc. Amer. 24 
[1952], 230. 

[2] A report on the First ICA Congress is given in J. acoust. 
Soc. Amer. 26 [1954], 110 and in Acustica 3 [1953], 390. 
The full proceedings, including all of the papers 
presented, are reprinted in Acustica 4 [1954], no. 1. 


EXPERIMENTELLE UNTERSUCHUNGEN 
AN SCHWIMMENDEN ESTRICHEN MIT SCHALLBRUCKEN 


von M. Heck 


Mitteilung aus dem akustischen Laboratorium von Prof. CREMER, Miinchen 


Zusammenfassung 


Nach einer Theorie von L. CREMER sind folgende GréBen fiir die Verbesserung der Trittschall- 
dimmung von schwimmenden Estrichen mit punktformigen Briicken maSgebend: Die Eigen- 
frequenz des Systems aus Estrichmasse und federnder Schicht, der Widerstand des Estrichs und 
der Rohdecke, die Ausbreitungsfunktion im Estrich und die Ubertragungsfunktion der Briicke. 
Es werden Verfahren angegeben, diese GréBen zu messen. Hierbei ergab sich der Absolutwert 
des Widerstandes und die Ausbreitungsfunktion einer 2 m2 Asphaltplatte in guter Ubereinstim- 
mung mit dem theoretischen Wert fiir eine unendlich groBe Platte. Die Rohdecke bestand aus 
einer Kappendecke mit lose dariiber vergossenem Zementglattstrich; deshalb unterschied sich 
der Widerstand sehr von dem fiir eine homogene Platte berechneten. Es wurden schwimmende 
Estriche mit kiinstlichen Briicken gebaut. Die gemessene Verbesserung stimmte gut mit der nach 
der obigen Theorie vorausberechneten iiberein, wenn man den tatsachlich gemessenen Wider- 
stand der Rohdecke einsetzt. 


Summary 


According to a theory of L. CREMER the following factors are important for the insulation 
against air-borne sound of floating floors with sound-bridges: the natural frequency of the system 
floating floor and resilient layer, the resistance of the floor and covering, the propagation function 
in the floor and the transmission function of the supports. Methods of measuring these factors are 
given and, as example of the consonance between the theory (for an infinite floor) and practice, 
the measured values of resistance and propagation function in a 2 m? asphalt floor are quoted. 
The measured improvement in floating floors constructed with special supports also agrees with 
theory when the resistance of the covering is included. 


Sommaire 


D’aprés une théorie de L. Cremer, les quantités importantes en ce qui concerne l’amélioration, 
au point de vue amortissement du bruit des pas, des carrelages flottants avec ponts acoustiques 
pontiformes, sont les suivantes: la fréquence propre du systéme composé de la masse du carrelage 
et de la couche élastique, la résistance du carrelage et du gros ceuyre du plafond, la fonction propa- 
gation dans le carrelage et la fonction propagation des ponts acoustiques. On indique des procédés 
de mesure de ces quantités. On a trouvé alors que la valeur absolue de la résistance et la fonction 
propagation d’une plaque d’asphalte ayant 2 m? sont en bon accord avec les valeurs théoriques 
relatives 4 une plaque indéfinie. Le gros ceuvre du plancher se composait d’une chape, sur laquelle 
avait été coulé un enduit lisse de ciment n’adhérant pas; c’est pourquoi la résistance était trés 
différente de celle qui a été calculée pour une plaque homogéne. On a construit des carrelages 
flottants, avec ponts acoustiques artificiels. L’amélioration mesurée cadre bien avec celle calculée 
a partir de la théorie ci-dessus, si on y introduit la résistance effectivement mesurée du gros 
ceuvre du plafond. 


1. Theorie 
Fiir die Verbesserung AL der Trittschalldam- 


bestimmten Ubertragungsfaktor, der folgender- 
mafen definiert ist [1]: | 


mung durch schwimmende Estriche mit punkt- n= 2o 1 F, 2 

formigen Briicken hat L. CremER [1] folgenden fing, bat wet (e+ 2) (2) 
Ausdruck abgeleitet: 

1 mit vy = Schnelle des Estrichs am Briickenort, 

AL = 10 log (fif)* + IZ, hea NP (1) wenn keine Briicke vorhanden ware, 0) = 


Schnelle der Briicke an dem dem Estrich zuge- 


Dabei bedeuten f die jeweilige Frequenz, f, die 
Eigenfrequenz des Systems Estrichmasse gegen 
Federwirkung der Zwischenschicht, Z, den 
Biegewelleneingangswiderstand des Estrichs, Z, 
den Eingangswiderstand der Rohdecke, /7 die 
Ausbreitungsfunktion fiir Biegewellen im Estrich 
und N den im wesentlichen durch die Briicke 


wandten Ende, v, = Schnelle der Briicke an dem 
der Rohdecke zugewandten Ende und F, = von 
der Briicke auf den Estrich zuriickwirkende 
Kraft. 

Fiir einen zu Longitudinalschwingungen ange- 
regten Stab mit der Dichte 0,, der Flache A, und 
der Elastizitat E,, dessen Lange h, klein ist gegen 
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die Wellenlange, wird in der obengenannten 


Arbeit folgender Ausdruck angegeben: 


wes Z, . {A101 Z, 
N= (1471) +i(G% +a )om @) 
Fiir den weiteren Spezialfall einer sehr weichen 
Briicke kann man in Gl. (2) das zweite Glied ver- 


nachlassigen, da in diesem Falle die riickwirkende 
Kraft sehr klein ist. Es gilt also: 


N ® 0,/vp. (4) 


Diese Annaherung erweist sich insofern als 
sinnvoll, als sie einen haufig auftretenden Fall 
wiedergibt, und weil dieser Quotient, den wir im 


folgenden mit (4a) 


bezeichnen wollen, nur mit seinem Absolutbetrag 
in Gl. (1) eingeht und dieser relativ einfach zu 
messen ist. Wollte man N genau messen, so 
miBte man auch die Phasenbeziehungen zwi- 
schen v, und F, erfassen, was meftechnisch 
einige Schwierigkeiten bereitet. 

Fiir die in den meisten Fallen verwendete An- 
regung tiber eine Masse m (Hammermasse) ist 
in Gl. (1) deren Widerstand zu beriicksichtigen. 
Hierzu ist Z) in Gl. (1) zu ersetzen durch Z, + 
iwm. Der zweite Summand im Nenner ist also 
zu multiplizieren mit Z>/(Zj + @?m?). Dasselbe 
gilt aber auch fiir den ersten Summanden, der 
zu dem iiber die elastische Schicht tretenden 
Energieanteil gehért [2]. Wir erhalten somit ein 


additives Glied und aus Gl. (1) wird: 


1 
+ 
(filf)* + |Z211/Z. N/? 
ZR + w? m? 
Zi 
Mit dieser Formel werden die in dieser Arbeit 
angefihrten MefBergebnisse verglichen werden. 


Doch dazu ist es erst nétig, die GréBen f,, Z,, Zo, 
IT und N zu kennen. 


Ay, = V/V, 


AL = 10 log 


(5) 


+ 10 log 


2. MeBobjekte 


Die verwendete Rohdecke hatte eine Flache 
von 10m? und eine Dicke von ca. 25 cm. Sie 
bestand aus einer Kappendecke, die nachtraglich 
mit einem Zementglattstrich von ca. 3 cm ver- 
sehen wurde. Leider war die Verbindung zwischen 
Kappendecke und Glattstrich sehr lose und auch 
ziemlich ortsabhangig. Da aber in Gl. (5) nur der 
Eingangswiderstand der Rohdecke eingeht, fiihrt 
dieser Mangel zu keiner Verfalschung der Meb- 
ergebnisse, wenn man die Briicke nur an den 
Stellen der Rohdecke aufsetzt, deren Eingangs- 
widerstand bekannt ist. 
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Als Estrich wurde eine Asphaltplatte von 2 m? 
Flache und einer Dicke von hy, = 2,1 em ver- 
wendet. Durch Messungen an einem Probestab 
aus einem moéglichst gleichartigen Material wur- 
den nach dem von R. Scumipt entwickelten Ver- 


fahren [3] folgende Materialkonstanten be- 
stimmt: 
Dichte o, = 2.0 g/cm", 
Longitudinalwellen- 


geschwindigkeit cro = 
Verlustfaktor 7 


2,2 - 10° cm/s, 
0,1 (bei 6°C). 

Der Zwischenraum zwischen Estrich und Roh- 
decke wurde mit Glaswolle ausgefillt, d. h. mit 
einem Material, dessen Anteil an der resultieren- 
den Federung gegeniiber dem der eingeschlosse- 
nen Luft nahezu vernachlassigt werden kann [4]. 

Man kann also f,; berechnen nach der Formel: 


fr = 000//Gd, (6) 
wobei G das Flachengewicht des Estrichs in 
kg/m? und d den Abstand Estrich—Rohdecke in 


cm bedeuten. 


I 


3. Messung der Eingangswiderstande 


Zur Messung der Eingangswiderstande diente 
ein K6rperschallsender, der aus einem Perma- 
nentmagneten der Masse mm = 10500 g und 
einer Schwingspule auf einem Aluminiumkopf 
der Masse m = 580 g bestand. Um den Frequenz- 
gang der Apparatur aufzunehmen, wurde der 
Magnet und der Schwingkopf auf eine sehr weiche 
Unterlage gestellt. 

FlieBt ein Strom I durch die Schwingspule, 
dann wirkt auf den Aluminiumkopf und auf den 
Magneten eine entgegengesetzt gleiche Kraft. 


F=aI=mbs=mm bm (7) 


(bs = Beschleunigung der Schwingspule, bu = 
Beschleunigung des Magneten). 

In Abb. 1 sind die mit dem Kérperschall- 
abtaster gemessenen Werte von 20 log (bm/b,) 
und 20 log(bs/b,) (6) = willkiirlicher Bezugs- 
wert) aufgetragen. Wie man sieht, ist der Ab- 
stand der beiden Kurven bis ca. 800 Hz konstant 


40 =, 
dB 
Abb. 1. oa 5 ae 
Gemessene Kor- 
perschallpegel 20 


bei Anregung 
mit konstantem 


Strom; bs auf 1 

dem Schwingkopf z a 

gemessen, by, auf tee | 

dem Magneten 0 100 200 400 9800 1600 3200Hz 
gemessen. ——> Frequenz 
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und zwar ist entsprechend Gl. (7) 
20 log ee: 20 ae a) 
bu m 


Uber 800 Hz ist der Abstand nicht mehr gleich, 
weil hier bereits Eigenschwingungen des Magne- 
ten auftreten; das heiBt also, das der Magnet 
nicht mehr als starrer Kérper anzusehen ist. DaB 
gleicher Strom nicht zu gleichen Kérperschall- 
pegeln fiihrt, diirfte auf einen Frequenzgang des 
Beschleunigungsabtasters zuriickzufiihren sein. 
Fiir diese Annahme spricht auch, daB man, halt 
man nach Abb.2 ein Kondensatormikrophon 
unmittelbar tiber die Schwingspule, bis 3000 Hz 
nahezu eine konstante Druckanzeige erhalt 
(Abb. 3). Dieser Druck aber diirfte bei der vor- 
liegenden Anordnung zum iiberwiegenden Teil 
von den Beschleunigungen der Luft in dem 
Zwischenraum herriihren, die radial gerichtet 
sind und radial zum Rande wachsen. Der Ab- 
stand betrug hierbei etwa nur 3 mm, ist also 
sicher klein zur Wellenlange. Aus dem gleichen 
Grunde dirften auch am Ringspalt, dessen 


Abb. 2. Anordnung zur Messung des Schalldrucks iiber 
dem Schwingkopf. 


~ 
o 
| 


a 
wo 


—-~= Druckpegel 


ra) 


| 
700 200 400 800 1600 3200He 


——= Frequenz 


Abb. 3. Schalldruck iiber dem Schwingkopf bei Anregung 


mit konstantem Strom. 
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Durchmesser 3,9 cm betrug, die Tragheitsanteile 
die Strahlung bei weitem iiberwiegen. Es ist da- 
her anzunehmen, da der gemessene Druck un- 
mittelbar ein Maf fiir die Beschleunigung der 
Schwingspule ist, und die Feststellung, daB dieser 
Druck konstant ist, l4Bt auch auf ein konstantes 
x und ein Verhalten der Schwingspule als trage 
Masse schlieBen. Zumindest diirften Anderungen 
von «, die an sich durch Anderungen der Form 
des Magnetfeldes méglich sind, von so geringer 
Bedeutung sein, da ihre Beriicksichtigung fiir 
bauakustische Messungen nicht lohnen wiirde. 
Da der Eingangswiderstand definiert ist als 


Z = Flvp=io F/bp (9) 
(F = erregende Kraft, vp = Geschwindigkeit, 
bp = Beschleunigung der Platte), und man bei 


einem iiber alle Frequenzen konstanten Strom 
die Kraft als konstant ansetzen kann, erhalt 


man: (10) 

Die Bestimmung der Konstanten erfolgt aus 
Gl. (7): (11) 
(b},;= Beschleunigung des Magneten unter 800 Hz). 
Also wird aus Gl. (10): 


Z=iwmmbi,/bp. 


Z = io const/bp. 


const = mm by, 


(12) 

Zur Messung ist es also nur nétig, bei kon- 
stantem Strom den Kérperschallabtaster einmal 
auf die Platte und bei einer Frequenz unter 800 
Hz auf den Magneten aufzusetzen und die erhal- 
tenen Werte entsprechend dem in Abb. 1 aufge- 
tragenen Frequenzgang zu korrigieren. Es ist 
dabei nicht notwendig, den Absolutbetrag der 
Beschleunigung und der Konstante « zu kennen. 

Weit einfacher wird die Messung unter 800 Hz, 
wenn die Gleichung 


F = mm bu (13) 
noch voll giiltig ist. Setzt man Gl. (13) in Gl. (9) 
ein, so erhalt man: 
, , bu 
Z=iwmm ret 

Man braucht also nur die Pegelunterschiede 
auf dem Magneten und auf der Platte zu messen 
und erhalt daraus direkt den gesuchten Wert. 
Im Prinzip vergleicht man den Widerstand der 
Masse des Magneten mit dem Widerstand der 
Decke. 

Da es nicht méglich war, die Phasenbeziehun- 
gen von bp und by zu messen, wurden in Abb. 4 
nur die Absolutbetrage der Widerstande einge- 
zeichnet. Messung a stellt den Eingangswider- 
stand der oben beschriebenen Asphaltplatte dar. 
Gestrichelt ist der theoretische Wert fiir eine 
unendlich groBe Platte eingetragen, der sich aus 
dem Widerstand der Platte [2] 


(14) 
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Abb. 4. 

(a) Eingangs- 
widerstand einer 
Asphaltplatte 
(gestrichelt: 
theoretische 
Werte); 

(b) Eingangs- 
widerstand einer 


Rohdecke. 
10° 
100 200 400 800 ‘600 3200Hz 
——> Frequenz 
Zi 2,9 .€Lo 0, he (15) 
und dem Widerstand iwm der erregenden 


Schwingspule zusammensetzt: 
Z| = 2,3? eh.ophs +o? m®. (16) 
Wie man sieht, ist die Ubereinstimmung sehr 
gut, wenn man die oben angegebenen Werte fir 
CLo, Qg und hy, einsetzt, besonders wenn man be- 
ricksichtigt, daB die Platte nur 2 m? gro} war. 
Aber wegen der hohen inneren Dampfung von 
Asphalt diirfte zumindest am Anregungsort die 
direkte Welle die reflektierten weit iiberwiegen. 
Die Streuungen in Abhangigkeit vom Anregungs- 
ort waren sehr klein, so da man annehmen kann, 
daB die angegebenen Werte fiir die ganze Platte 
giltig sind. Messung b stellt den Widerstand der 
Rohdecke dar. Hier treten groBe Unterschiede von 
den nach Gl. (16) berechneten Werten auf; dies 
diirfte aber mit der obenerwahnten ungeniigen- 
den Verbindung zwischen Kappendecke und 
Zementestrich zusammenhaingen. Nimmt man 
an, daf nur der Glattstrich erregt wurde, so er- 
halt man einen den Messungen entsprechenden 
Wert von Z,. Ob diese Annahme zu Recht be- 
steht, ist fiir das weitere belanglos, denn es wurde 
immer mit den tatsichlich gemessenen Werten 
von Z, gerechnet. 


4, Ausbreitungsfunktion im Estrich 


Zur Messung der Ausbreitungsfunktion // 
wurde der Estrich mit dem Kérperschallsender 
angeregt und einmal direkt am Anregungsort und 
dann in verschiedenen Entfernungen die Be- 
schleunigung gemessen. Die so erhaltenen Werte 
sind in Abb. 5 aufgetragen. Gestrichelt wurden 
die Werte eingezeichnet, die man fiir eine unend- 
lich groBe Platte zu erwarten hatte [3]. 


(kr) =H? (kr) —H(—ikr) — (17) 
H = Hankelsche Funktion, k = komplexe Wel- 
lenzahl, r = Abstand vom Anregungsort. Wie 
man sieht, kann man noch fiir Abstande bis zu 


70 em die Asphaltplatte angenahert als unend- 
lich ausgedehnt betrachten. 
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Abb. 5. 
Ausbreitungs- 
funktion der 
Asphaltplatte 
(gestrichelt : 
theoretische 
Werte). 


1700 200 400 800 
——* Frequenz 


1600 3200Hz 


5. Ubertragungsfaktor der Briicke 


Fiir stabférmige Briicken interessieren die 
GréBen A,, 0,, h, und E,. Die ersten drei sind 
leicht zu bestimmen, zur Messung von FE, wurde 
folgendes Verfahren angewandt (Abb. 6). Auf 
den Schwingkopf der Masse m wird die zu unter- 
suchende Briicke und darauf die Masse M ge- 
klebt. Der Sender wird mit reinen Ténen (Kraft 
F) angeregt und auf M wird die Beschleunigung 
gemessen. Dann ist: 


Abb. 6. Anordnung zur Messung des Elastizitatsmoduls 


der Briicke. 
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F,—F,=iomn,, (18) 
Sale, 
Us == 05 ee ee 19 
1 2 E, A, 2 ( ) 


F, = F,+iw0,A,h,v,, F,=iw Mv, (20), (21) 
(F,, v, sind Kraft und Schnelle am Anfang, F,, v, 
sind Kraft und Schnelle am Ende der Briicke). 
Aus Gl. (18) und (21) erhalt man: 
elec et a Po « (22) 
iw(m+ M+A,h,0,—@?m Mh,/E, A,) 
Daraus erhalt man fiir die Resonanzfrequenz: 
E, A, 
mMh, 


(m+ M + A, hy, 0;) (23) 


oRmMhy 
A,(m+ M + A, hy Q,) ; 


Fiir die spater verwendete Korkbriicke ergab 
sich z. B. E, = 1,3 - 108 dyn/cm?. 

Die fiir schallweiche Briicken mafgebende 
GréBe A,. wurde folgendermafen bestimmt 
(Abb. 7): A,, hangt von den Briickeneigenschaf- 
ten und von Z,, nicht aber von Z, ab. Man setzt 
also die Briicke auf die Rohdecke, auf die spater 
der Estrich gebracht werden soll, auf und regt 
die Briicke mit dem Kérperschallsender an. Dann 
mift man einmal die Beschleunigung am oberen 
Ende der Briicke b4 und dann am unteren Ende 


br. Da (25) 


ist, hat man damit bereits den gesuchten Wert. 
In Abb. 8 sind einige so gemessene Werte |/,,| 
aufgetragen. Die untersuchten Briicken waren: 

1. ein Gummi-Element (Abb. 9), 

2. eine sogenannte Biigelbriicke aus Eisen 

(Abb. 9), 

3. eine Korkbriicke. 

Beim Gummi-Element ist deutlich die erste 
Eigenresonanz zu erkennen, die nach anderen 
Messungen etwas iiber 800 Hz liegen miiBte. Bei 
der Bigelbriicke liegen die Verhaltnisse kompli- 


On 


und fir EH: £,= (24) 


Ay, = v4/Vz = by /bz 


FS 
i 
ae 


Abb. 7. Anordnung zur Messung des Faktors A,,. 
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dB 
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— > 20 log Ay) 
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Abb. 8. 
MeBwerte des 
Faktors Aj, fiir 
einige Briicken; 


(1) Gummibriicke, 
(2) Bigelbriicke, 100 200 400 800 1600 3200Hz 
(3) Korkbriicke. ——> Frequenz 


Abb. 9. Gummi-Element und Biigelbriicke. 


zierter. Auf jeden Fall ist A, viel kleiner, als 
nach dem statisch bestimmten Elastizitatsmodul 
des Biigels zu erwarten ware. Wahrscheinlich 
erfolgt die Ubertragung durch Biegeschwingun- 
gen des Eisenbiigels. 


6. Messungen an schwimmenden Estrichen mit 
kiinstlichen Briicken 


An der Stelle, an der die Briicke angebracht 
werden sollte, wurden erst der Eingangswider- 
stand und der Trittschallpegel der Rohdecke be- 
stimmt. Dann wurden die Briicke und der Estrich 
aufgesetzt und der Zwischenraum mit Glaswolle 
angefillt, so daB der Abstand Estrich-Rohdecke 
konstant war. Da es sich gezeigt hatte, daB die 
Verbindung der Briicke nach beiden Seiten hin 
von groBem EinfluB ist, wurde die Briicke oben 
und unten angeklebt. Zur Anregung diente teils 
das normale Hammerwerk, teils der Kérper- 
schallsender. 

Bei dieser Gelegenheit wurde auch versucht, 
bei der Rohdeckenmessung einen Absolutver- 
gleich der beiden Anregungsarten durchzufiihren. 
Dieser Vergleich beruht darauf, daB, wachst der 
Widerstand der Decke weniger stark mit der 
Frequenz an als der Widerstand der auffallen- 
den Masse, dies physikalisch damit gleichbedeu- 
tend ist, dafS die Masse mit der Decke fest ver- 
bunden ist und auf die Masse eine dem Stof ent- 
sprechende Kraft wirkt. Im vorliegenden Fall ist 
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100 200 400 800 1600 3200Hz 


——> Frequenz 


Abb. 10. Trittschallpegel bei Anregung mit dem Hammer- 
werk und mit dem KG6rperschallsender. MeBwerte 
auf gleiche Kraft und Bandbreite korrigiert. 

(a) Hammerwerk, 
(b) K6rperschallsender, 
(c) Differenz. 


der Widerstand der Decke einigermaBen konstant 
und die Masse des Schwingkopfes und des Fall- 
hammers sind annahernd gleich; obige Voraus- 
setzungen sind also sicher erfiillt. 

Abb. 10 zeigt die gemessenen Trittschallpegel 
bei Anregung mit dem Hammerwerk, wenn man 
die Werte auf eine Bandbreite von 75 Hz korri- 
giert. AuBerdem zeigt sie auch die Trittschall- 
pegel, die man erhalt, wenn man die Decke mit 
dem Kérperschallerreger mit einer sinusférmigen 
Wechselkraft von 1,7-107’ dyn anregt. Diesen 
Wert erhalt man, wenn man das vom Hammer- 
werk erzeugte Kraftspektrum berechnet und die 
Energien innerhalb des oben angegebenen Berei- 
ches summiert unter der Annahme, da die ein- 
zelnen Teilténe inkoharent sind. In Abb. 10 
wurde auch der Unterschied der beiden An- 
regungsarten angegeben, fiir dessen Frequenz- 
gang Th. LANGE [5] einen ahnlichen Verlauf er- 
halt. Wenn man beriicksichtigt, da es sich hier 
um eine Absolutmessung unter ziemlich undefi- 
nierten Voraussetzungen (Lagerreibung des Ham- 
merwerkes und dergleichen) handelt, ist die er- 
haltene mittlere Abweichung von 2,8 dB als 
gering zu bezeichnen. Fiir Relativmessungen ist 
dieser Unterschied natiirlich zu groB. Es wurden 
daher fiir die Verbesserung nur die Werte gleicher 
Anregungsart herangezogen. 

Abb. 11 zeigt die Verbesserung bei Vorhanden- 
sein einer Eisenbriicke von h, = 1,1 cm und 4,= 
16 cm? in verschiedenen Abstanden von der An- 
regungsstelle. Gestrichelt ist die nach Gl. (5) ge- 
rechnete Verbesserung eingetragen, wenn man 
die gemessenen Werte fiir Z), Z,, N und die ge- 
rechneten Werte von // zugrunde legt. Der fir 
r = 0 ersichtliche Anstieg von AL fiir hohe Fre- 
quenzen kommt nicht von einer Federwirkung 
der Briicke, sondern vom Eingangswiderstand 
der vorgeschalteten Masse m (zweites Glied in 
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100 200 400 800 1600 3200Hz 
——> Frequenz 


Abb. 11. Verbesserung bei Vorhandensein einer Eisen- 
briicke in verschiedenen Abstéanden vom An- 
regungsort (gestrichelt = theoretische Werte). 


Gl. (5)). Die verwendete Briicke wirkt also als 
vollkommen starre Verbindung. Da man aus 
diesen Messungen schlieBen kann, daB die Ent- 
fernungsabhangigkeit in der Formel richtig er- 
faBt wurde, wurden die folgenden Messungen nur 
mehr direkt an der Anregungsstelle vorgenom- 
men. 

In Abb. 12 ist die gemessene Verbesserung als 
ausgezogene Kurve eingetragen. Als Briicke 
wurde ein stabférmiges Korkstiick mit folgenden 
Daten verwendet: A, = 20cm’, h, = 1,9 cm, 
0, = 0,2 g/cm’, E, = 1,3 - 108 dyn/em?. Die ge- 
strichelte Kurve stellt die Verbesserung dar, die 
man erhalt, wenn man die obigen Daten in Gl. (3) 


30 


| 


| 
100 200 400 800 


——= Frequenz 


Abb. 12. Verbesserung bei 
briicke; 
———w— nach Gl. (5) und (3) berechnete Werte 

fiir eine kurze, stabférmige Briicke; 

nach Gl. (5) und (4a) berechnete 

Werte unter Zugrundelegung der ge- 

messenen Werte von Aj); 

—+- +: — gemessene Werte. 


1600 3200Hz 


Kork- 


Vorhandensein einer 
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und (5) einsetzt. AuBerdem ist in die Abb. 12 
punktiert die Verbesserung eingetragen, die sich 
ergibt, wenn man in Gl. (5) und (4a) die in Abb. 8 


aufgetragenen Werte von A,, fiir Kork einsetzt. 


Wie man sieht, sind die drei Kurven einiger- . 


mafen gleich, nur die aus A,, berechnete punk- 
tierte Kurve liegt erwartungsgem4B etwas nied- 
riger, weil in diesem Fall die auf den Estrich 
zuriickwirkende Kraft vernachlassigt wird. 

Wie man auBerdem aus Abb. 12 ersieht, zeigt 
die vorliegende Bricke bei tiefen Frequenzen 
keine Wirkung, weil hier die Schallenergie haupt- 
sachlich durch die elastische Schicht tritt, auf der 
der Estrich schwimmt. Man erhalt die normale 
Verbesserung: 


AL = 40 log (f/f,): (26) 


Erst bei héheren Frequenzen macht sich der 
EinfluB der Briicke bemerkbar; die Verbesserung 
steigt nicht mehr so stark an, weil die Briicke als 


Federung wirkt. Wie man aus Gl. (3) und (5) 


schlieBen kann, gilt fiir héhere Frequenzen: 
AL = 20 logaf 


(a ist eine bestimmte Konstante). 

Abb. 13 und 14 zeigen die Verbesserungen 
beim Vorhandensein des Gummi-Elementes bzw. 
der Bigelbriicke. Die ausgezogenen Kurven be- 
deuten wieder die gemessenen und die gestrichel- 
ten Kurven die nach Gl. (3) und (5) berechneten 
Werte, wenn man fiir das Gummi-Element den 
nach Gl. (24) berechneten Elastizitatsmodul EF, = 


(27) 


100 200 400 800 


——> Frequenz 


1600 3200Hz 


Abb. 13. Verbesserung bei Vorhandensein einer Gummi- 
briicke (Bedeutung der Kurven siehe Abb. 12). 


40 
dB 


100 200 400 800 
——= Frequenz 


Abb. 14. Verbesserung bei Vorhandensein einer Biigel- 
briicke (Bedeutung der Kurven siehe Abb. 12). 


1600 3200 Hz 


M. HECKL: SCHWIMMENDE ESTRICHE MIT SCHALLBRUCKEN 


ACUSTICA 
Vol. 5 (1955) 


1,1 - 107 dyn/cm? einsetzt. Bei der Biigelbriicke 
war diese Messung nicht méglich, deshalb wurde 
E, aus der statischen Durchbiegung bestimmt. 
In beiden Fallen unterscheiden sich die gemesse- 
nen und so gerechneten Werte von AL ganz er- 
heblich. Das hat verschiedene Ursachen: Wie 
schon erwahnt, ist fiir das Gummi-Element die 
Voraussetzung, daB die Lange h, klein ist gegen 
die Wellenlange, nicht mehr erfillt, weil schon 
bei 800 Hz die erste Kigenresonanz auftritt. DaB 
in diesem Frequenzgebiet die Verbesserung nicht 
noch weiter absinkt, diirfte auf die hohe innere 
Dampfung von Gummi zuriickzufihren sein. Bei 
der Biigelbriicke hangt die Verbesserung iiber- 
haupt nicht mehr mit dem statisch bestimmten 
Elastizitatsmodul zusammen, denn Ahnlich wie 
ein dinner Blechstreifen, hat der Eisenbiigel 
schon bei tiefen Frequenzen (90 Hz) Eigenreso- 
nanzen und leitet deshalb die Schwingungen in 
Form von Biegewellen sehr gut von der Estrich- 
platte auf die Rohdecke. 

Bei der Messung von A,, und der daraus be- 
stimmten Verbesserung (punktierte Kurve) wer- 
den diese Einfliisse mit erfaBt. Daher unterschei- 
den sich in diesem Fall die gerechneten und ge- 
messenen Werte nur wenig. Die Vernachlassigung 
der Riickwirkung der Briicke bedeutet bei diesen 
sehr weichen Bricken keinen groBen Unterschied, 
denn je weicher die Briicke ist, um so kleiner ist 
auch die auf den Estrich zuriickwirkende Kraft. 
Fiir den allgemeinen Frequenzgang der Verbesse- 
rung gilt auch bei diesen beiden Briicken das 
oben Gesagte: Wenig EinfluB der Briicke bei 
tiefen Frequenzen, groBe Verschlechterung bei 
den héheren. 

Die vorstehenden Untersuchungen wurden im 
Auftrage und mit Mitteln des Bundes-Ministe- 
riums fir Wohnungsbau durchgefihrt. 

Herrn Professor Dr. L. CREMER bin ich fiir die 
Anregung zur vorliegenden Arbeit und seine zahl- 
reichen wertvollen Ratschlage bei ihrer Durch- 
fiihrung zu groBem Dank verpflichtet. 


(Eingegangen am 25. August 1954.) 


Schrifttum 


[1] Cremer, L., Acustica 4 [1954], 273. 

[2] Cremer, L., Acustica 3 [1953], 317. 

[3] Scumipt, R., Ing. Arch. 5 [1934], 352. 

[4] GosrLx, K. und Bacu, W., Mitteilung des Instituts fiir 
Technische Physik, Stuttgart. 

[5] Lancr, Tu., Acustica 3 [1953], 161. 

Siehe auch: 

[6] Exxinc, W., Acustica 4 [1954], 396. 

[7] Cremer, H. und L., Frequenz 2 [1948], 61. 

[8] GéseLe, K., Uber die Wirkungen von schwimmenden 
Estrichen (unveréffentlicht). 

[9] Meyer, E., PARKIN, P.H., Operst, H. und Purkis, H. J., 
Acustica 1 [1951], 18. 


ZUR THEORIE ALLGEMEINER, 


SCHWACH GEKOPPELTER AKUSTISCHER RESONATOREN 


von E.. CHRISTIAN 
Institut fiir Theoretische Physik der Universitat Graz 


Zusammenfassung 


Im Sinne der Anwendung von Niherungsverfahren auf allgemein nicht lésbare Randwert- 
probleme der Elastizitatstheorie wird eine Stérungsrechnung erster Ordnung auf ein schwach 
gekoppeltes verlustfreies akustisches Resonatorensystem angewandt. Im allgemeinsten Fall 


handelt es sich um n Resonatoren, die durch (3) Koppelglieder (die teils wegfallen kénnen) ener- 


getisch verbunden sind. AuBer der Frequenzaufspaltung infolge der Verstimmung durch die 
Koppelraiume wird sie durch Volumendeformation berechnet. Bei beliebiger Gestalt der Begren- 
zungsflachen des Systems wird fiir diese ein fester bzw. freier Rand vorgeschrieben. Schwiichere 
Randbedingungen und einschrinkende Bedingungen beziiglich des Schwingungstyps werden im 
Rahmen von gesonderten Betrachtungen beriicksichtigt. Neben schwacher Kopplung der Resona- 
toren durch die Koppelriume und einer gemeinsamen Eigenfrequenz der Einzelresonatoren im 
Ausgangszustand setzt die Rechnung keine wesentlichen Einschrinkungen voraus. Die Koppel- 
frequenzen ergeben sich als Lésungen einer Sakulardeterminante, deren Elemente durch Energie- 
terme gebildet sind. Als Spezialfall wird ein System mit schubspannungsfreiem Medium behandelt. 
Der funktionelle Zusammenhang der Kopplungsparameter eines Zweiresonatorensystems mit den 
Resonanzerscheinungen wird klargelegt und abgeleitet. SchlieBlich werden einige praktische 
Beispiele durchgerechnet. 


Summary 

A perturbation calculation of the first order is made on a weakly coupled loss-free system of 
resonators. Beside the spread of frequency due to mistuning by coupling, that due to volume 
deformation is also reckoned. For given boundary conditions rigid or flexible walls are prescribed. 
More and less restricted conditions with respect to the mode of vibration are considered. In respect 
of weak coupling through the cavity and general values of the frequencies of the individual 
resonators, the calculation makes no special restriction. The coupling frequencies emerge as 
solutions of a secular determinant of which the elements are made up of energy terms. As special 
case a medium free of shear is considered. The functional relationship of the coupling parameter of 
a double resonator is derived. Finally some practical examples are given. 


Sommaire 

Pour appliquer un procédé d’approximation 4 la résolution de problémes de la théorie de 
Vélasticité relatifs aux valeurs aux parois et qui ne sont généralement pas solubles, on effectue 
un calcul de perturbations du premier ordre sur un systéme de résonateurs acoustiques sans pertes 
et faiblement couplés. Le cas le plus général est celui de n résonateurs reliés énergétiquement par 


(2) éléments de couplage (dont une partie peut manquer). On calcule la multiplication de la fré- 


quence produite par le désaccord dfi aux cavités de couplage et celle provenant de la déformation 
en volume. On attribue au systéme, dont les surfaces limites ont une forme quelconque, une paroi 
rigide ou une paroi libre. On traite des cas particuliers ot les conditions aux limites sont moins 
sévéres et moins limitatives quant au mode d’oscillation. Le calcul ne comporte aucune restriction 
importante, en dehors de la condition de couplage lache des résonateurs par les cavités de couplage 
et de fréquence propre commune 8 tous les résonateurs. Cas particulier d’un systéme avec milieu 
exempt de contraintes de cisaillement. Relation fonctionnelle entre les paramétres de couplage 
d’un systéme a deux résonateurs et les phénoménes de résonance. Quelques exemples pratiques. 


1. Einleitung 


Das vorliegende Problem der Elastizitats- 
theorie ist zu den stehenden Schwingungen, im 
Gegensatz zu den in unendlichen Medien sich aus- 
breitenden Wellen, einzuordnen. 

Zur rechnerischen Erfassung beliebiger Schwin- 
gungen eines allgemeinen akustischen Resona- 
torensystems ist eine strenge Integration der 
elastischen Gleichungen unter Beriicksichtigung 


des Randverhaltens der Feldvektoren an den Be- 
grenzungsflachen des aus Resonatoren und Kop- 
pelraumen bestehenden Gesamtraumes durch- 
zufiihren. Diese Integration ist aber nur fiir geo- 
metrisch besonders einfache Hiillflachen (Réhre, 
diinner Stab, Kreisring, Platte usw.) méglich. 
Die klassischen Arbeiten iiber gekoppelte 
Schwingungen von n elastischen Resonatoren 
gehen im wesentlichen auf Lord Ray .eicu [1] 
guriick, Auch die Arbeit von WIEN [2] iiber die 
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Riickwirkung eines Resonators ist in diesem Zu- 
sammenhang erwahnenswert. In einer zusam- 
menfassenden Darstellung tiber gekoppelte Re- 
sonatoren faBt BackHaus [3] n mit einander ge- 


koppelte Einzelsysteme als ein System mit n- 


Freiheitsgraden auf. 

Von GrieBE und BLEcHsScHMIDT [4] wurde eine 
spezielle Kopplungstheorie fiir ein zylindrisches 
Einresonatorensystem entwickelt. Von der Kop- 
pelschwingungsgleichung ausgehend, wurde von 
den Verfassern eine Kopplung der radialen und 
axialen Eigenschwingung angenommen. Sie konn- 
ten ihre Rechnung allerdings nur fir Hohl- 
zylinder geringer Wandstarke experimentell be- 
statigen, wahrend sich fiir Vollzylinder keine 
Ubereinstimmung mit der Theorie ergab. 

Eine formale Analogie zur Theorie der elektri- 
schen Doppelleitung stellt die akustische Lei- 
tungstheorie dar. Sie ist fiir Widerstandsbe- 
stimmungen, im Zusammenhang mit Anpassungs- 
fragen und Kopplungsproblemen, besonders aber 
zur Erfassung der MeBgréBen des Schallfeldes 
von besonderer Bedeutung. In der Folge wird 
daher kurz auf die MeBgréBen der akustischen 
Leitung eingegangen, um den spater sich als not- 
wendig erweisenden Ubergang zu den Feld- 
gréBen der Theorie elastischer Schwingungen zu 
finden. 

Wird ein akustisches Leitungselement (Stab, 
Trichter, Hohlleiter usw.) durch eine Eingangs- 
und Ausgangsflache charakterisiert, so stellt es 
den Typ des akustischen Vierpols dar. Der Ab- 
schlu8 des Vierpols ist in zwei Extremfallen be- 
sonders bedeutungsvoll: 

a) der schallharte Abschlu8 (Schallschnelle ver- 
schwindet), 

b) der schallweiche AbschluB (Schalldruck ver- 
schwindet). 

Diesen Extremfallen entspricht der akustische 

Eingangswiderstand bei Leerlauf bzw. bei Kurz- 

schluB. 

Die den elektrischen LeitungsgréBen Wechsel- 
spannung und Wechselstrom entsprechenden 
GréBen sind Schallwechseldruck und Schall- 
schnelle. Der Betrag der Impedanz ist gleich dem 
Verhaltnis von Druck und Schnelleamplitude. 
Die Phase der Impedanz ist die Phasendifferenz 
von Druck und Schnelle. Ist der Druck an der 
Begrenzungsflache konstant (Rohrleiter), so laBt 
sich die Impedanz zweckmaBiger aus Druck und 
SchallfluB definieren. So ist die Schnelle in einem 
luftgefiillten Rohrleiter als konstant anzusehen 
(ebene Welle) und sinkt nur in der Reibungszone 
der Wand auf Null. In diesem Fall ist der FluB 
dem Produkt Querschnitt mal Schnelle im Rohr- 
zentrum gleichzusetzen. Im Falle einer unendlich 
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langen Rohrleitung stellt das Verhaltnis Druck 
zu Schnelle den Schallwellenwiderstand dar. 


2. Der homogene akustische Leiter 


Da sich der Hauptteil der Arbeit auf verlust- 
freie Schwingungen eines elastischen Systems be- 
zieht, sollen auch in diesem leitungstheoretischen 
Teil nur verlustfreie Energieleitungen beriick- 
sichtigt werden. Die homogene akustische Leitung 
stellt die Analogie zu dem symmetrischen elektri- 
schen ‘Vierpol dar. Sie wird als verlustfrei be- 
zeichnet bei verschwindendem Widerstandsanteil 
R und bei Verschwinden der Ableitung G. Es 
verbleiben demnach nur die rein imaginadren 
Widerstandsanteile. Diese Voraussetzungen tref- 
fen etwa bei einem mit reibungs- und schub- 
spannungsfreiem Medium ausgefiillten Metallrohr 
von konstantem Querschnitt und geniigender 
Wandstarke zu. Die SchwingungsgréBen der aku- 
stischen Leitung sind der Schalldruck Q und der 
Schallflu8 J. Als weitere wichtige MeBgréBe wird 
der Wellenwiderstand Z eingefihrt. Fir die 
Schwingungsgréfen gilt die Beziehung: 


6/dt = jo (2,1) 
Fir den Wellenwiderstand Z gilt: 
Z=VRiRc, Z—) Ruka) 


wobei mit Rw der Massenwiderstand und mit Rr 
der elastische Widerstand bezeichnet wird. Das 
FortpflanzungsmaB g wird dargestellt durch: 


g=a+jb=VRijRc, g=\VRu/Re. (2,3) 


Wegen der vorausgesetzten Verlustfreiheit ver- 
schwindet das DampfungsmaB a. Mit 6 wird das 
die Phasenverschiebung bestimmende Winkel- 
maf} bezeichnet. Die SchwingungsgréBen ergeben 
sich fiir eine beliebige Stelle x der Leitung nach 
der akustischen Vierpoltheorie als: 


Qs = Qa coshg x —Ja Z sinhg x, (2,4) 


a 


(2,5) 


Jx = Ja coshg x —-> sinhg x. 
Die Indizes a weisen darauf hin, daB es sich bei 
diesen GréBen um die Anfangsbedingungen han- 
delt, wobei fiir die Stelle mit dem belastenden 
Widerstand Q./J. die Koordinate x = 0 gesetzt 
wird. Fiir die dampfungsfreie Leitung ergibt sich 
weiter: 

Qs = Q. cosh x —j Ja sinbx, (2,4’) 


(2,5’) 


Jz=Jacosbx —j Qe sinbx. 


Z 
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Die Transformationsgleichung der Impedanz 
lautet fiir eine beliebige Stelle x der Leitung: 


pee i Ztgbt og) 
Jz  1—(R,./Z)-tgbx 
An Hand dieser Beziehung laBt sich durch Hinzu- 


schalten eines definierten Vierpols ein Widerstand 


R. in den Widerstand R, transformieren. Mit 
Hilfe eines Vierpols kann die Messung von Im- 
pedanzen fiir verschiedene Schichtdicken bei vor- 
gegebener Frequenz oder bei vorgegebener 
Schichtdicke fiir verschiedene Frequenzen durch- 
gefiihrt werden. 


3. Grundlagen, Rand- und Zusatzbedingungen 


Die Entwicklung einer allgemeinen Theorie 
schwach gekoppelter akustischer Resonatoren 
setzt einige Ergebnisse der Elastizitatstheorie 
voraus. Effekte, bei denen die Schallabsorption 
eine Rolle spielt, bleiben unberiicksichtigt. Den 
Ausgangspunkt bildet die fiir sehr kleine Ver- 
rickungen eines homogenen und isotropen Kér- 
pers giiltige elastische Grundgleichung: 

be 100 
Os > 
0a a — div P, (0 1 ) (3,1) 
001 


Verwendete Bezeichnungen: 


0 = Dichte im Ruhezustand, 

A, = Lamésche Konstanten, 

w= anf = Kreisfrequenz, 

§ = (u,,U,,u,) = Vektor der elastischen 
Verschiebung, 

Ss = Deformationstensor, 

P = Spannungstensor, 

div P = Vektordivergenz P, 

div P = skalare Divergenz P, 

G = Gradiententensor. 

Das Medium selbst wird durch die beiden 


Laméschen Konstanten « und / charakterisiert. 
Der Spannungstensor hat die Gestalt: 


P=2yuS+AdivsE. (3,2) 


Die Vektordivergenz des Spannungstensors P 
lautet: 
div P= wAS+(u+A)graddivs. (3,3) 


Bei Beriicksichtigung der Synchronitatsbedin- 
gung 


= 3 f(t) Fei 


(3,4) 
fiir den Verschiebungsvektor, gewinnt die elasti- 
sche Grundgleichung (3,1) die zweckmaBige Ge- 


stalt: —o w®3 = wAS+(u+A)graddivs. (3,5) 
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An Randbedingungen, die sich sowohl auf Re- 
sonatoren als auch auf Koppelelemente der Kon- 
figuration beziehen kénnen, werden vorgeschrie- 
ben: 


=) 3,6 

a ; fiir den festen Rand, ) 

(s, 72) = 0 (3,6a) 
Sayey eT 
Cow fir den freien Rand. (3,7) 
(Pn, n)=0 (3,7a) 
) 


Die schwacheren Forderungen (3,6a) und (3,7a 
sind nur bei einer Spezialisierung auf schub- 
spannungsfreie Medien zuladssig. Eine Hinzu- 
nahme von Zusatzbedingungen kommt einer 
Einschrankung der allgemeinen Lésungsmannig- 
faltigkeit gleich. Oft sind es allerdings derartige 
Zusatzbedingungen, die die Lésung als physi- 
kalisch besonders sinnvoll erscheinen lassen. 

Als Zusatzbedingungen sollen beriicksichtigt 
werden: 


a) Einschrinkende Bedingungen beziiglich des 
Feldverlaufs in den Resonatoren 


Derartige Zusatzbedingungen kommen einer 
Spezialisierung des Schwingungstyps gleich. Wird: 
1. divs = 0 gesetzt, so erhalt die Grundgleichung 
die einfache Gestalt: 


0: 63/2 = wAs. (3,8) 


Die Lésungen der Gleichung (3,8) bezeichnet man 
wegen des Fehlens jeglicher Kompression als 
Scherungswellen. 

2. rots=0 gesetzt, so ist der Verschiebungs- 
vektor als Gradient eines skalaren Potentials dar- 
oe rots = —gradq. (3,9) 
Die Grundgleichung laBt sich in eine Wellen- 
gleichung 


0: &p/dt? = (24+ A) Ag (3,10) 


des Potentials p iiberfiihren. Sie vereinfacht sich 
mit “= 0 fiir ideale Fliissigkeiten und Gase noch 
weiter auf: 


0: @p/a2=AAg. (3,10a) 


b) Stetigkeitsbedingungen fiir geschichtete Medien 
Die zweite Gruppe der Zusatzbedingungen sind 
Stetigkeitsbedingungen fiir die Grenzflachen zwi- 
schen zwei verschiedenen Medien. Ausfiihrliche 
Betrachtungen tiber Reflexion und Brechung an 
ebenen Grenzflachen wurden von ScHocu [5] an- 
gestellt. Die den drei Fallen: 
1. zwei schubspannungsfreie Medien, 
2. festes und schubspannungsfreies Medium und 
3. zwei feste Medien grenzen aneinander 
entsprechenden Stetigkeitsbedingungen lauten: 
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(P, n,n) =(P.,n) ae 
4 G,, 2) = (pn) Cae 
(Pasn) = (Ea) eee 
a G,, 7) = (Sy, 7) { (3, ) 
3. ate (3,11e) 


Wendet man das Prinzip von actio und reactio 
an [6], so kann gezeigt werden, dass die Steitig- 
keitsbedingungen dem Energieprinzip geniigen. 


4. Allgemeine storungstheoretische Entwicklung 


Zur Berechnung der Koppelfrequenzen eines 
allgemeinen, schwach gekoppelten n-Resonato- 
rensystems wird eine Stérungsrechnung erster 
Art angewendet. 

Die durch (3) Koppelglieder verbundenen n 
Resonatoren sollen im entkoppelten Zustand 
mindestens eine gemeinsame Frequenz besitzen. 
Die Einzelresonatoren seien R, mit v = 1,2,---,n, 
wahrend die Koppelglieder mit R,;, r= 1, 2, ---,n; 
s=1, 2, ---, (n —1), bezeichnet werden. Das Me- 
dium des Systems sei reibungsfrei, fest und iso- 
trop. Die Ausgangsfrequenz wm, wird von den 
Eigenfrequenzen der Koppelraume genigend 
weit enifernt angenommen. Im Ausgangszustand 
sind die Resonatoren entkoppelt. Ein stetiger 
Frequenziibergang von w,) zu zu @, wird sodann 
durch das Offnen von in den Koppelflachen K,. 
angebrachten Blenden bewerkstelligt. Die Fre- 
quenzen @, sind demnach Eigenwerte der Feld- 
verteilung des Gesamtraumes vom Volumen V5 
im Endzustand. Es werden folgende Bezeichnun- 
gen verwendet: 


V, = Volumen des v-ten Einzelresonators, 
V,; = Volumen des rs-ten Koppelraumes, 
F, = Flache des v-ten Resonators, 

F,; = Flache des rs-ten Koppelraumes, 


Ry = R+ YR; = Gesamtraum, bestehend 
el C{r, s) 
aus den Resonatoren und 


den Koppelraumen R,., 


= 


Vs =¥ V+ ¥ Vs = Gesamtvolumen des Sy- 
ies Cirs| stems, 

Fy =) F,+ Y¥ F,. = Oberflache des Gesamt- 
poy OA systems. 


Das Koppelglied R,, stellt das Verbindungsstiick 
zwischen R, und R, dar. Die Eigenwertgleichung 
des Systems ist sodann von der allgemeinen 


Gestalt: . 
f(A, Ms ©, Pyy +s Pn) = 0. (4,1) 


Neben dem Eigenwertspektrum enthilt sie impli- 
zit die Konstanten 4 und yw sowie die Gestalts- 
parameter pj mit 1= 1, 2,---,n. Es wird weiter 
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von Vorteil sein, auch fiir den ungestérten Zu- 
stand einen geringen energetischen Kontakt 
innerhalb des Gesamtsystems vorauszusetzen. 
Diese Voraussetzungen geniigen bereits, um das 


‘ Problem erfolgreich stérungstheoretisch behan- 


deln zu kénnen. Im Ausgangszustand handelt es 
sich um eine n-fache Entartung, denn es fallen 
die n Eigenfunktionen $, v= 1, 2, ---,n, jeweils 
nur mit der Feldverteilung des v-ten Resonators 
zusammen, wahrend sie in allen anderen den 
Betrag Null haben. 

Die n Eigenfunktionen gehen im Sinne der 
Stérungsrechnung nicht stetig in jene der ge- 
stérten Feldverteilung tiber. Dies gilt aber fiir 
eine den $s, zugeordnete lineare Schar 

S45) enue dey (4,2) 
os 
sobald man die Koeffizienten a geeignet be- 
stimmt. Sowohl den ungestérten als auch den 
gestorten Feldverlauf driicken wir durch die in 
beiden Fallen giiltige elastische Grundgleichung 
aus. Es gilt daher: 


—w 0% = div Py, (4,3) 
—w* 9 3, = div Pi. (4,4) 

Fir die variierte FeldgréBe wird gesetzt: 
Sk = Sk + OS: (4,5) 


Weiter kann man im Sinne einer Stérungsrech- 
nung erster Ordnung setzen: 
O=W)+60, w=ai+2a,dw, (4,6) 
P;, = P,.+ oP. k= 1, 2s a (4,7) 
Die Indizes k beziehen sich auf die dazugehérige 
Eigenfunktion §,. Aus den Gleichungen (4,3) und 


(4,4) erhalt man nach Einfiihrung der variierten 
GréBen: 


div 6 Pi + w2 0 8% = —2 wy) dm 0%. 


mit 


(4,8) 


Multipliziert man die letzte Gleichung der Reihe 
nach skalar mit den Divergenzen der n Span- 
nungstensoren, so folgt daraus: 


div P\? div Px + w2 0 0% div Pi? 


= —2 bw si div P™, 4,9 
0 


Es seien 37 (v=1,2,---,) die Eigenfunk- 
tionen eines energetisch abgeschlossenen Systems, 
jeweils gebildet aus dem Resonator R, und den 
unmittelbar angrenzenden Koppelgliedern Ros 
(s=1, 2, ---, n—1). Die den Eigenfunktionen ${7) 
entsprechenden Tensoren werden mit P(?), 6, 
S\%) bezeichnet. 

Fiir das Teilsystem R, z ys Re lautet die 
Grundgleichung: 
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rela 
Abgeschlossenes System 


(n=1) 
R, +: > ea ° 
s=1 
eingefiihrt zur Definition der 
EKigenfunktionen $5) 32). 


= (2) 


div P? + w?293!7 =0, (4,10) 


wiahrend fir die variierten Gréfen gilt: 


39 = 3,468, w2=w2 +2azdo. 


Damit erhalten wir aus Gl. (4,10) das Gleichungs- 
system: 


div PR? = —@2 03% —2.0,6wa.5). (4,11) 


Multipliziert man Gl. (4,11) jeweils mit div6 P; 
skalar, so erhalt man im Zusammenhang mit 
Gl. (4,10) bei Vernachlassigung von Gliedern die 
klein sind von héherer Ordnung, das aus n* 
Gleichungen bestehende System: 
wo (6%, div PX”) = —2 w,dwo0(%, div P!”’) + 

+ wo (8, divd Px). (4,12) 


Durch Verwendung der Identitat: 


div Ps = (3, div P)+ PG (4,13) 
erhalt man weiter: 
wo (div P\? 63, — P\” 6 Gi) 
= 2 dwo (&, div div P”) + 
+ we o(divé P.3%—5P.C). (4,14) 


Nach Integration tiber das Gesamtvolumen unter 
AusschluB der Koppelraume, erhalt man bei 
gleichzeitiger Anwendung des GauBschen Satzes: 


2.0600 [ (%, div P\”) dt + 
Vn— 2 Vea 


Cir, s} 


+a [( bP, G2— PMG)a 


Da ZX Vrs 


C{r,s) 


(4,15) 
= wo fle P37, n— P63, n) df + 


Fp 


ote T PORK, 
Cir, s| F, 


n— P 6%, n) df. 


Die letzte Gleichung 1aBt sich unter Beriicksichti- 
gung der Randbedingungen (3,6), bzw. (3,7) nach 
einfacher Rechnung weiter vereinfachen auf: 
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os [(6 P,3?) ’ n)— (P P35, n)| df 
C{r, s} EF 
Xa) (Z) 
S93 (5, div P’) d (4,16) 


In dieser Gleichung sollen nun durch weitere Um- 
formungen saimtliche gestérten Feldgréfen eli- 
miniert werden. So erhalt man bei Beriicksichti- 
gung der giltigen Beziehungen: 


= (2) 


+(Z 
a, 36 (Fab 


Davee Sie ere oe =3$ +05, 
3 SO =z Ou 

sowie der Randbedingungen (3,6), baw. (3,7) nach 
Vernachlassigung von kleinen Gliedern héherer 
Ordnung bereits die einfache Form: 


— $( (P? 3 3, n) df 


— 


C{r, s} Fes 
i oe [ Gx. div P,) dr. (4,17) 
pare VE Yn 
{r,s} 


Nach erneuter Anwendung des Gaubschen Satzes 
als auch der Identitat (4,13) erhalt man das 
Gleichungssystem in der gewiinschten Form: 


fo div ee POG. | dz 


(4,18) 


Fiir die einzelnen Integrale, die Energieterme 
darstellen, verwenden wir die Abkiirzungen: 


v= if (div P., 5,) dz, (4,19) 
hae >» Vrs 
Cir, 5} 

dae if [(div P, 3?) + POG] dz. (4,20) 
Vrs 

Wegen _ ye Ok p 3° 

p=1 
und C= grad a = Xkp grads)” 


erhalt man bei Verwendung der Abkiirzungen 


(4,19) und (4,20) das Gleichungssystem: 
{((aiv P®, 3?) + POG?) dr = ph aap AN’). 
Vrs ‘(4,21) 


Schreibt man die linke Seite von Gleichung (4,18) 
in der Form: 


[ (div P”, 3?) + PMG? dr 


(4,22) 
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und beriicksichtigt, da fiir die Terme im einzel- 
nen gilt: 


Bure 1. (r,s) ==; p) Le, ( 
P 2. (153) =, p) = ACE ==Acn ; 
f; 3. r==0v Gaye 
a red 4. r=v Spi we 4 


so laBt sich das erhaltene Gleichungssystem wei- 
ter vereinfachen. Die Abb. 2 veranschaulicht im 
besonderen die Terme der Gruppe (a, 2). 


R, Rrs Rs 
2(2)_-e(2) = 
SG?) Heise a 
' 


| 


Abb. 2. Energieterme eines 
Koppelgliedes der 
Gruppe (a, 2). 


Fiir die Energieterme W,, der Resonatoren 
selbst findet man: 


mit v-+p wird W,p>=0 
mit v=p wird W,,=—W. 
= 2 | (div P., 3.) dr. 


V>— » Vrs 
Cir, s) 


Durch die Verwendung der Abkiirzungen (4,19) 
gewinnt die rechte Seite von Gl. (4,18) die 
Gestalt: 


2 if (5, div P,) dt == a Xkp Wp = Xkv Wr . 
p=l1 


V>— DZ Vrs 
C{r, s} 


Desgleichen schreibt man statt Gl. (4,22): 


2d Xkp S Ae 


p=1 Cir, s} 


n n 
/ , v,s 
= 3 Xkp Aup + Oke ye Ae ) . 
p=. s=1 


Damit erhalt das grundlegende Gleichungssystem 
die iibersichtliche Form: 


yy onp Arp + Xke eee 2) 9@ W..| = 0 (4,23) 
p=1 s=1 Wo 

(US 1)2, s,s bl 25:-, n)- 
Dieses lineare Gleichungssystem besitzt als 


Lésungen der als Sakulargleichung bekannten 
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Determinante ebenfalls symmetrisch. Setzt man 
den Rang der Determinante mit (n—1) voraus, 
so folgt daraus die Existenz von n reellen Lésun- 
gen. Es sind dies die erstrebten Aufspaltungs- 
frequenzen Sao 

Die Sakulargleichung ist allgemein anwendbar, 
da sie als RechengréBen ausschlieBlich Schallfeld- 
eroBen der ungestérten Konfiguration enthalt'. 
Sie wurde unter Beriicksichtigung sowohl der 
festen als auch der freien Berandung abgeleitet. 
Von den Termen sind die Aww fiir die Kopplung 
verantwortlich, wahrend die A’ die Riickwir- 
kung der energetisch. abgeschlossenen Koppel- 
glieder auf die einzelnen Resonatoren bestimmen. 

Zusammenhange elementarer Natur, wie sie 
fiir die Sakulargleichung Giltigkeit haben, sind 
vor allem fiir die Falle der Praxis von Bedeutung. 
Es ist offensichtlich, daB sich die geometrische 
Anordnung der Resonatoren und ihrer Koppel- 
glieder in der Struktur der Determinante wieder- 
spiegelt. Es ist dies der Fall fiir Resonatoren oder 
Koppelglieder gleicher Beschaffenheit, Serien- 


schaltungen usw. 


5. Stérungsrechnung des skalaren Problems 


Die allgemeine schubspannungsfreie Materie 
gestattet nicht mehr die vereinfachende Annahme 
der Inkompressibilitat. Die Dichte wird an jeder 
Stelle des Raumes durch den dort wirkenden 
Druck in eindeutiger Weise bestimmt. Der 
Spannungstensor entartet zu einem Skalar, wo- 
durch die elastische Grundgleichung die wesent- 
lich vereinfachte Form 


o 07s /ét? = grad P 
annimmt. Mit 


(5,1) 
P=) divs E 
gilt fiir die Vektordivergenz P: 
div P = 7 grad (div3) = grad P. 
Nach den im Abschnitt 3 gebrachten Ergeb- 


nissen gelangt man nach Einfiihrung des skalaren 
Potentials zu der Wellengleichung 


Systemdeterminante Agt+k?o=0, worin k?=o9w?/d 
Meal woe Ansghek tone 7 gesetzt ist. Das hier vorliegende 
s ee a Problem ist auf ein skalares zuriick- 
An, ~ An! —60Wx ... Ar ang gefiihrt. Durch sinngemaBe Wieder- 
: : ‘ ne holung der stérungstheoretischen 
Ate Asatte. eee > A® 0 EBS, Pr ee (4,24) 1 Fiir schwach gekoppelte verlustfreie 


die gesuchten Aufspaltungsfrequenzen. Wegen 
der Symmetrie der Terme A,,—Ap ist die 


elektromagnetische Hohlraumresonatoren 

kommen E. LepinecG und P. URBAN 
(Acta Phys. Austriaca 2 [1948], 198) zu einer ahnlich ge- 
bauten Sikulardeterminante. 
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Uberlegungen soll die spezielle Form der Sakular- 
determinante fiir schubspannungsfreie Medien 
gewonnen werden. 

Die Eigenfunktionen, fiir die ein stetiger Uber- 
gang in den Randbedingungen bestehen soll, 
werden hier ebenfalls als Linearaggregate 


Pr = ys Xkv Pov 
v=) 
mit k=1, 2,---,n dargestellt. 
Fiir die variierten GréfBen gilt 
Pk = Pk+OGe, W=W_t+dw, w?=w2+ 20,dw. 


Die Wellengleichungen, sowohl fiir den unge- 
stérten als auch den gestérten Zustand angesetzt, 


(5,2) 
(5,3) 


Agr = —(9.0§/A) vx, 
AGr = —(Q.0?/A) Gr, 


bilden den Ausgangspunkt. 

Wegen der sinnvollen Wiederholung der ein- 
zelnen Uberlegungen des vorangegangenen Ab- 
schnittes, die rechnerisch gestaltet auBerdem eine 
wesentlich einfachere Form besitzen, wird hier 
auf die Teilschritte verzichtet. 

Die der Gl. (3,18) entsprechende System- 
gleichung lautet fiir das skalare Problem: 


Sf (gradg! grad gf” + 9 Ag) ax 


Cir, s} Vrs 
= |v Aq dt, 


PS V7 
2 C{r, s) 


(5,4) 


wobei die Indizes der Eigenfunktionen sinnge- 
maf iibernommen wurden. Dieses Ergebnis stellt 
daher einen Spezialfall des allgemeinen vekto- 
riellen Problems dar. 


6. Das vektorielle Problem bei Verstimmung 
durch Voluniendeformation 


Es wird im folgenden das Problem erweitert. 
Fir den Fall einer Verstimmung durch Deforma- 
tion wird die Frequenzgleichung des Resonatoren- 
systems abgeleitet. 

Auf Grund der Tatsache, dafs ein schwach 
gekoppeltes akustisches Resonatorensystem un- 
ter bestimmten Voraussetzungen einen Vierpol 
darstellt, 1aBt sich die abgeleitete Theorie zur Er- 
fassung der Filterwirkung eines derartigen Vier- 
pols verwenden. Zur Bestimmung der Durchla- 
gebiete einer schwach gekoppelten Resonatoren- 
konfiguration sei eine gegenseitige Verstimmung 
méglich. Diese wird im wesentlichen durch eine 
Volumendeformation erreicht. Die Eigenfre- 
quenzen sind gleichzeitig Grenzfrequenzen der 
DurchlaBbereiche. 
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Zur Beibehaltung der vorangegangenen Me- 
thode wird eine gemeinsame Eigenfrequenz der 
Einzelresonatoren im entkoppelten Zustand an- 
genommen, die Kopplung des Systems durch einen 
simultanen Vorgang bewirkt und schlieBlich die 
Verstimmung durch eine Deformation erzwun- 
gen. Die Kopplungsanderung vom Nullwert zum 
Endwert wird durch das Offnen von in den Kop- 
pelflachen angebrachten Blenden realisiert. 

Die Rechnung lauft iibereinstimmend mit Ab- 
schnitt 4 bis zur Beziehung (4,14). 

Zur Auswertung der Randbedingungen (3,6) baw. 
(3,7) wird sodann iiber das Volumen des Endzu- 
standes unter AusschluB der Koppelraume R;. 
integriert. Das durch die Deformation der Reso- 
natoren verdrangte Volumen wird mit AV, be- 
zeichnet, wahrend F's die Hiillflache des Gesamt- 
systems im Endzustand darstellt. Durch Inte- 
gration tiber das Volumen 

Vice AV 
Cr, s} 7a 

und gleichzeitige Anwendung des Gaufschen 
Satzes erhalt man aus Gl. (4,14) das System: 


ye (3, divP) dr + 
Oy vas SAV, 
at iGArzsh rw 


a f (6 P, 6G — P\5Gx) dr 


n 
V>— 3) Vrs— JX AVr 
C{r, s} r=1 


=D $[(6 P30, n) — (P53, n)] af+ 


C{r, s} Fe, 


(6,1) 


ag f [(6 Pk 30, AWS (Pi 0% F n)| df. 


Fy 


Zur Gewinnung einer Ubersicht des Integrations- 


bereiches der Gl. (6,1) dient die Abb. 3. 


Abb. 3. Integrationsbereich der Gl. (6,1). 
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Der Integrand in der zweiten Zeile von Gl. (6,1) 
verschwindet exakt, wahrend in den der ersten 
Zeile 


>) >( +(Z) 


= 5x OS; =e + OSz 


eingefihrt wird. 


Fiir das Integral in der vierten Zeile von 


Gl. (6,1) wird gesetzt: 
J (Fs) = [ (8 Ped, n)— (PP 9%, n)] df. (6,3) 


Tae 
Fir die Flache Fy gilt die Randbedingung 
(S;, n) = 0. Mit F, wird die aus Fy und Fy zu- 
sammengesetzte Flache bezeichnet. In dem Inte- 


gral a 
J(F;) =J(Fs) + J(Fs) = J(Fs) 


verschwindet J(F'y) wegen der Randbedingung. 

Substituiert man Gl. (6,2) in (6,3), so bekommt 
man bei Vernachlassigung kleiner GréBen héherer 
Ordnung nach einfacher Rechnung fiir das Inte- 
gral in der vierten Zeile von Gl. (6,1): 


(6,4) 


JF) => f [(div P,, %) + P»GiJdr. (6,3’) 


r= 


AV, 

Der Summenausdruck in der dritten Zeile von 
Gl. (6,1) 148t sich unter Beriicksichtigung der 
Randbedingung (3,6) durch einfache Rechnung 


umformen auf: 


Side Ps siaar 


Cirs) yy 


= Yl, div PX?) + PL? G1?) ac. 


C{r, s) Vrs 


(6,5) 


Damit sind samtliche gestér- 


ten FeldgréBen aus Gl. (6,1) 


eliminiert und sie erhalt die An, 

Gestalt: : 
20 =: nis 
“2 (div P., 54) de 


{(@?, div P®) + PGP]dr+ (6,6) 


[(&. div P,) + Py Gx] dr. 


asta 

Zu dem gleichen Ergebnis gelangt man fiir die 
Randbedingung (3,7). Die Rechnung verlauft 
analog, lediglich an Stelle von Gl. (6,2) ist die 
entsprechende Beziehung der Spannungstensoren 
zu verwerten. 

Zur Gewinnung einer besseren Ubersicht wer- 
den Abkiirzungen eingefiihrt. Bei stetigem Uber- 
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rs s é6 
yi AY ay Peete | 58 
s Wo 
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gang vom Anfangs- zum Endzustand wird das 
passend gewahlte Linearaggregat der Eigen- 
funktionen 
Sh ¥" Oke So, [Fi Wp oaeih 
v=1 
stetig in den Endzustand ubergehen. Durch die 
Einfihrung der analogen Beziehung 
3 = Y ons, kk 
p=! 


== 1,2, soon 


in Gl. (6,6) lassen sich die Energieterme durch 
die vorgegebenen Feldvektoren ausdriicken. Die 
Abkiirzungen (4,19) und (4,20) werden erganzt 
durch: 

B= Hs [(div P.,3,)+ P»Gp] dt. (6,7) 
AV 

Es gilt: 

ra } rT) ve) 

Bw r=p=. 


Damit geht die Beziehung (6,6) in das lineare 
System: 


n ae 6 
yy okp ‘Alp + OXkv (ras ) ae Ba Ws} aa 0 
p= M9 


(v= 1,2, - ne kk = ho 2 (6,8) 


uber. Den hinzugekommenen Termen B,, kommt 
die Funktion der Verstimmung des Systems zu. 
Mit allen Termen B,.=0 geht das Problem in 
das in Abschnitt 4 behandelte iiber. 

Das in den «xp lineare und homogene System 
(6,8) ist durch das Verschwinden der System- 


determinante 


Aj Nea e Ain 
Ago pO Ow Aon 
: = 2 
A A v4 B _ bo W 
n2 - “inn nn Wo nn (6,9) 


in nichttrivialer Weise lésbar. Die n Eigenfre- 
quenzen kénnen unter Umstinden zusammen- 
fallen, doch sind sie wegen der Symmetrie- 
bedingung Aix, = Axi stets reell. 


7. Kopplungsparameter eines 
Zweiresonatorensystems 


Die Methode der Verwendung elektrischer 
Ersatzschaltbilder hat sich in der akustischen 
Leitungstechnik bewahrt. Auch der Bau akusti- 
scher Siebketten hat den elektrischen Analogie- 
fall zum Vorbild. Durch Gl. (6,9) lat sich neben 
den méglichen Koppelschwingungen auch die 
Filterwirkung eines derartigen Systems angeben. 
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Sie enthalt bereits die Grenzpunkte der Durch- 
laBbereiche. 

Doch auch die Kopplungsparameter lassen 
sich daraus allgemein ableiten, was in diesem 
Abschnitt an Hand eines Zweiresonatorensystems 
gezeigt wird. Dieses Ergebnis ist vor allem fiir 
die auf dem Gebiete der Ultraakustik im Vorder- 
grund stehenden Fragen der Anpassung, Kopp- 
lung und Verstimmung im Zusammenhang mit 
dem Energietransport durch akustische Leitun- 
gen [7] von Bedeutung. 

Es wird nun der funktionelle Zusammenhang 
der Kopplungsparameter mit den Energietermen 
fiir ein symmetrisches Zweiresonatorensystem 
abgeleitet. Der allgemeine Charakter der Ablei- 
‘tung beruht auf der Verwendung der Feld- 
vektoren. 

Als Kopplungsparameter werden eingefiihrt: 
a) der Kopplungsgrad k, 

b) die Verstimmung des Einzelresonators p, 
c) der Verstimmungsgrad des Gesamtsystems v. 

Die Kopplungsparameter sind dimensionslos 
und erweisen sich fiir die durch die energetischen 
Wechselwirkungen der Einzelraume bewirkte 
Frequenzaufspaltung verantwortlich. Die Saku- 
lardeterminante lautet fiir ein symmetrisches 
Zweiresonatorensystem : 


Ow 
Ae eee Wy, Ars 
Wo 


An, Agy) + Bo2— ae Wo» 
0 

Wird nur ein Resonator verstimmt, was prinzi- 
piell keine Einschrankung bedeutet, es sei R,, so 
kann man B,,=0 setzen. Die Energieterme von 
Gl. (7,1) setzt man durch die physikalisch begriin- 
deten Kopplungsparameter an. 

Wenn W,, = W,,. = W die Gesamtenergie eines 
Einzelresonators ist, so laBt sich der gemischte 
Term des Koppelraumes mit dem Verstimmungs- 
parameter proportional setzen: 


Air” = Aly = pW. 


Weiter wird der aus den FeldgréBen eines ein- 
zelnen Nachbarraumes sich zusammensetzende 
Energieanteil des Koppelraumes, mit dem Kopp- 
lungsgrad proportional gesetzt: 


46=4,,=kW 
SchlieBlich wird die Energie des Deformations- 
raumes mit dem Verstimmungsgrad proportional 
gesetzt: he Sat 
Den Zusammenhang zwischen Koppelfrequenzen 


und Systemparametern erhalt man durch Ein- 
setzen der letzteren in Gl. (7,1). 
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Die Lésungen der Gleichung zweiten Grades, 

(dw /w»)?—(2 p+ v)da/w)+ p?+ pu—k?=0, 
sind die Aufspaltungsfrequenzen 

(Swi /o9) = > (2p+utfo?+ 4k); i=1,2. 


Vernachlassigt man den Verstimmungsanteil des 
Koppelraumes, so bekommt man mit p=0 die 
Aufspaltungsfrequenzen 


dwilw) =~ (vt Vo? + 4k), i=1,2. 
2 


Den Inhalt dieser Gleichung veranschaulicht 


Abb. 4. 


Abb. 4. Aufspaltungsfrequenzen als Funktion 
der Verstimmung v, Parameter ist der 


Kopplungsgrad k. 


Der Kopplungsgrad k ist Scharparameter, fiir 
k=0 wird 6@/m, eine lineare Funktion von v. 

Auch der Fall ohne Verstimmung durch 
Volumendeformation ist von Interesse. Die Ver- 
stimmung wird nur durch die Koppelleitung be- 
wirkt. Die Koppelfrequenzen lauten in diesem 
Fall: 

(dw:i/wo) =ptk. 


Verschwindet zusatzlich k, so folgt (d@/@ )1,2 = p- 
Es existiert eine einzige Koppelfrequenz, Fiir 
p= 0 werden die Aufspaltungsfrequenzen gleich 
dem Kopplungsgrad. 

Wegen der allgemein gemachten Vorausset- 
zung 6w/m <1 wird hier ersichtlich, daB p und k 
ganz allgemein sehr kleine GréBen sein miissen. 
Sind die Koppelfrequenzen gegeben, indem sie 
zum Beispiel einer Messung zuganglich sind, so 
lassen sich die Parameter p und k mittels der 
symmetrischen Grundfunktionen ausdriicken als: 


p= (6@, + 6a) /2@o, k = (6@, + dm.) /2 m9. 
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8. Infinitesimale Deformation eines 
Einresonatorensystems 


Die Anwendungsméglichkeiten der Stérungs- 

rechnung auf ein akustisches Resonatorensystem 
sind vielseitig. Die Frequenzanderung eines in- 
finitesimal deformierten, schubspannungs- und 
reibungsfreien Einresonatorensystems kann vek- 
toriell (Abschnitt 4) oder skalar (Abschnitt 5) 
abgeleitet werden. Ist wy) ein nicht entarteter 
Eigenwert der Eigenfunktion $= 5%) = gradq@), 
so erfolgt von diesem Ausgangszustand ein Uber- 
gang zu einem benachbarten, durch Deformation 
gestérten Zustand. 

Die Rechnung ergibt fiir den festen und freien 
Rand dasselbe Ergebnis. Die vektorielle Stérungs- 
rechnung ergibt fiir die Frequenzanderung des 
Einresonatorensystems — es eribrigt sich auf 
die Teilschritte einzugehen, da es sich um eine 
starke Vereinfachung des allgemeinen Falles 
(Abschnitt 4) handelt — ausgedriickt durch das 
ungestérte Feld: 
mn [160 div Po) + PoSo] dt 


(8,1) 
Wo 2 | Poy ae 


Das skalare Rechenverfahren erbringt als gleich- 
wertiges Ergebnis, da es sich um ein schub- 
spannungsfreies Medium handeln soll: 


{(®o> div gradg, + gradqy, gradq,) dt 


6m av 


ayy genes Po, gradqyp) dt (8,1’ 


Anwendung 


Radialschwingung eines Kreiszylinders 

Es soll als Anwendungsbeispiel die Frequenz- 
anderung eines in der radialen Grundschwingung 
erregten, schubspannungsfreien Kreiszylinders 
nach Gl. (8,1) und (8,1’) berechnet werden. 

Fiir den Kreiszylinder sind die radi- dé 
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Von der Wellengleichung ausgehend, gelangt 


man zu der partikularen Lésung: 


p=Im = (=)",| cosm x Cos = Pa 


Mit m=0 und n=O erhalt man die radiale 
Grundschwingung g)=J,)(kr) und k=y,,/r 
= Y¥,/r, wobei unter yj, die erste Nullstelle der 
Ableitung der Besselschen Funktion nullter Ord- 
nung zu verstehen ist. Der Verschiebungsvektor 
hat fir die Lésung die Komponenten 


So{=koJi(kor), 0, 0} 


Es sind daher auch nur die Randbedingungen 
fiir den Mantel von Interesse. Fiir die Mantel- 
flache vor der Deformation (Radius =r) gilt die 
Randbedingung: 


OH o/ Or = —ky Jo (kor) = 0. 


Bei Verwendung von elementaren Funktional- 
gleichungen und Differentialformen fir Bessel- 
sche Funktionen lassen sich die Integrale der 
Gl. (8,1) einzeln auswerten. 

Man erhalt nach elementarer Rechnung: 


4 [Go div Py) + P, So] dt 


Wo 2 [ PS, dz 
we 


— R°[S5 ko R) + So (Ho R) So (ho R) —2 7° Sb (Kor)] 
= 27232 (kor) 


Nach Reihenentwicklung der Funktion 
f(R) = FR? [J5 (ho R) + Jo (ho R) Jz (ko R)] 


erhalt man weiter, bei Vernachlassigung kleiner 
GréBen héherer Ordnung, durch einfache Rech- 


nung: 


r2 Jz (kor) +2rdrJe(kor) —rdorJe(kor) —r? J (kor) 


alen Dehnungsschwingungen sowohl fir q@, 
feste isotrope Stoffe, als auch fiir schub- 
spannungsfreie Medien bekannt [8]. 
Wenngleich die Stérungsrechnung im allgemeinen 
auch erst bei exakt nicht lésbaren Problemen an 
Bedeutung gewinnt, so ist dieses Beispiel trotzdem 
besonders geeignet, da man eine zusatzliche Uber- 
prifungsméglichkeit an Hand der exakten Rech- 
nung hat. Der Einfachheit halber sei die Volumen- 
anderung durch die Anderung des Radius bewirkt. 

Verwendet man ein allgemeines Zylinder- 
koordinatensystem, so lassen sich die radialen 
Dehnungsschwingungen in bekannter Weise durch 
Besselsche Funktionen darstellen. 


7 Jo(kor) 


= Or/r. 


Als zweiter Schritt folgt die Berechnung der Fre- 
quenzanderung nach der skalaren Theorie (Gl. 
(8,1’)). Die Eigenfunktion der radialen Grund- 
schwingung lautet: g)= Jo(kor). Ihre Ableitun- 
gen sind gegeben durch: 


gradgy=S,, divgradg,= divs,= Ag, 
AGo= — kb Po = —k5 So (Kor). 


Durch elementare Auswertung der bestimmten 
Integrale erhalt man fiir die Frequenzanderung 
nach Gl. (8,1’): 
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te) col : 

bo Brioni mtgwl+ tga Ik (9,4) 
Wo j _—mtgoltgol, 

R? [Jo (ho R) + Jo (ko R) Ja (ko R)| —2 7? Jh (kor) 


272 J2 (kor) 


Als weitere Kontrolle wird die durch die exakte 
Rechnung erhaltene Eigenfunktion nach Taylor 
entwickelt. Bricht man nach dem zweiten Glied 
ab, so erhalt man bei Beriicksichtigung der Rand- 
bedingung: 


Ji (k R) =J, (kor) + Ok r Ji (kor) + dr ko Ji (kor), 
6k [ky = dw /w. = Or/r. 


9. Koppelschwingungen eines 


Zweiresonatorensystems 


Als weitere Anwendung der Theorie sollen die 
Koppelschwingungen zweier konzentrischer Lei- 
ter, die durch einen koaxialen Koppelleiter ener- 
getisch verbunden sind, berechnet werden. Das 
Gesamtsystem ist demnach eine dreifach zusam- 
menhangende, stiickweise homogene, akustische 
Leitung. Im klassischen Sinne handelt es sich 
hierbei um zwei schwach gekoppelte, einseitig 
offene Pfeifen. Die Abb. 5 veranschaulicht dieses 
Zweiresonatorensystem. 


x=0 ir 
0,=0 R,=R, Ry=R3 03=0 
R, Rie Ro 
ie ee (ko, I x 
q Z 


2 3 
le lo vat l3 i 
Abb. 5. Schnitt einer stiickweise homogenen Leitung. 


Nach Gl. (2,6) erhalt man fiir die Transformation 
der Impedanzen durch die Leitungslange: 


R,(1,) =—jZ,tgol,, (9,1) 
— R,—jZ.tewl 
Ry (lL, + 1) = : Lee 
1—j(R,/R,) tgo |, 
. (Z,/Z,) tg 1, + tga |, 
zee * 9,2 
14, 1—(Z,/Z,) tgol, tgol, oy 


Fiir die Koppelflache K,, 1a8t sich der trans- 


formierte Widerstand auch schreiben: 
R, (1, + 1s) =j Z,tgol,. (9,3) 


Die beiden letzten Gleichungen ergeben bereits 
die Eigenwertgleichung: 


a tye, — (AlZs) 80h + te0l, 


= . (9,4 
Z2 3 1—(Z,/Z,) tgol, tga l, (9,4) 


Fiir den Spezialfall der symmetrischen Anord- 
nung lautet die Eigenwertgleichung: 


wobei |, die Lange der Koppelleitung und Z;, den 
Wellenwiderstand der Koppelleitung bezeichnet. 
Weiter wurde in Gl. (9,4’) an Bezeichnungen ver- 
wendet: 
Vif: Z,.=25= 2: Z{Z.= ljk= im. 
Bei Anregung von R, in der axialen Grund- 
schwingung wird wegen w)y=22/A=a/l und 
®=W)+dqm die Frequenzengleichung bei Ver- 
nachlassigung kleiner GréBen weiter vereinfacht 
Zu: 


: 2 dw 1 
(de)? Se Lene = 0. (9,5) 
Die Lésungen 
61,2 = 2 (1 =f: : (9,6) 
“  Itg@ok \ ~~ cosagk y 


sind die gesuchten Aufspaltungsfrequenzen. 
Als zweiter Weg zur Auswertung dieses Bei- 
spiels soll nun der iiber die Sakulardeterminante 


beschritten werden. Sie lautet fiir ein Zwei- 
resonatorensystem: 
Ay? —s0Wn, Ais ~0. 07 
=0. (9, 
An, AS ) $6 Wo 
Wegen der Symmetriebedingung gilt: 
Ay? = AS”, Wu=We, Ap=An. 


Die SchwingungsgréBen an der Stelle x der Lei- 
tung ergeben sich nach der Leitungstheorie zu: 


2, 
(:) — Q.coswX; 


: ~ Oa .. 
(1) =—j Zi SIN W)X. (9,8) 
Wegen 


J, (2) =Jacoswgl=Ja, J, (l) =—j se SIN Wy lk 


Li 
ist weiter 
jZu 


. in ae 
sin Wolk 


0. = 


(9,9) 


Aus Gl. (9,8) und (9,9) erhalt man fiir die Lei- 


tungsgréBen des nach R, offenen Koppelraumes 


12° 


| Zed gies 
of = 1 Re eS, WX; 
sin @ kk 
J ee sin WX. (9,10) 
Sin Wo lk 


Ebenso erhalt man, durch Spiegelung und an- 
schlieBende Translation des Koordinatensystems, 
fiir das Teilsystem R,+ R,,: 
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(Z;) jZr Jo 
O° = intel COS (Ik —X), 
(Z,) Ie . ll 
J: ral gril sin@o(k—x). (9, ) 


An Hand dieser Beziehungen lassen sich nun die 
Integrale elementar auswerten: 


Awe Jel 


2 ° 


Wy =2 | (div P,,3,) dz= 
Vy 


70) 


Ay = J [(div PI®, 38) + PI? GI) de 
Vig 


ora Je 
2 Wy SiNWglk ” 


AG? = |[(Gv P, 30) + PO GM] ax 


12 
wean Je 


= CLE Oath 
WW 2 


Die Sakulargleichung (9,7) lautet demnach: 
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Fiihrt man wieder die Wellenwiderstande ein, so 
lauten die Aufspaltungsfrequenzen: 


Lik 1 1 
COLT 57; Tiga (1 sak) 


Damit ist volle Ubereinstimmung mit der nach 
der Leitungstheorie erhaltenen Lésung (9,6) er- 
reicht. 

Meinem hochverehrten Lehrer, Herrn Prof. Dr. 
P. UrBAN, sowie Herrn Prof. Dr. E. LEpINEGG 
bin ich fiir die Anleitung der Arbeit sowie das 
entgegengebrachte fordernde Interesse zu gréB- 


tem Dank verpflichtet. 
=n (Eingegangen am 4. Juni 1954.) 
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Bookreview 


E.Skuprzyk. Die Grundlagen der Akustik; 1st ed. Springer- 
Verlag, Wien, 1954. 1084 pages, 450 fig., 16 x 24 cm?, 
price DM 147,— or sfr. 150.—. 


The book gives in a concise form what scientists in 
acoustical research ought to know about theoretical 
acoustics. It is mathematical in the first instance, but over 
and again the author shows by practical examples how the 
theoretical results can be used advantageously for solving 
practical problems. 

The order in which the subjects are treated is rather 
unusual but, on closer examination, very well chosen. 
Fourier series, -integrals, Laplace transforms and transient 
phenomena are treated first, thus telling the reader that 
the knowledge of the steady-state behaviour of a system 
is sufficient for the treatment of all complex phenomena, 
but warning him at the same time to look out for pitfalls 
that he may not be aware of. 

The author proceeds then with the derivation and solu- 
tion of the differential equations of the sound field in rect- 
angular, cylindrical and spherical coordinates, exempli- 
fying their use in application to important problems such 
as reflexion, standing wave tubes, radiation, characteristic 
frequencies, diffraction, ete. Special chapters are devoted 
to electromechanical analogies, mechanical systems with 
distributed mass and compliance, coupled systems, trans- 
ducers, microphones, transmitters, piezoelectricity, speech 


and hearing, wave theory in room acoustics, room acoustics, 
sound absorption, lined ducts, dissipation of heat, under- 
water acoustics, etc. In the last chapter non-linear pheno- 
mena are dealt with. 

Relatively few references are given in the text, so that 
it is not always clear where the original publications are to 
be found. The book ends however with 126 pages of 
references, classified according to subjects in the same 
order as the text. This bibliography seems invaluable. 

Even when bearing in mind that many good books on 
acoustics have appeared in the last five years, I do not 
hesitate to state that the book fills a need. It is undoubtedly 
outstanding and on the whole astonishingly up-to-date. 
One might have doubts about the relative importance of 
subjects dealt with and the absence of certain topies, but 
it cannot be denied that each chapter gives so much 
valuable information that it is an excellent basis for 
further research. 

In view of the many outstanding qualities of the book 
it would be ungrateful to point to shortcomings and errors 
of minor importance that could be put right in a possible 
second edition. In spite of the rather high price I am 
convinced the book will find its way to all centres of 
acoustical research. 

As was to be expected in a book of a high scientific level, 
the MKS-system of units is used exclusively. 

C. W. Kosten 


DIE EINSTELLUNG VON TEILCHEN UND LOSUNGEN 
IM FELD STEHENDER ULTRASCHALLWELLEN 


von GERTRUD KECK 


Institut fiir Medizinische Physik der Universitat Wien, Vorstand: Prof. HAuER 


Zusammenfassung 


Die vorliegende Arbeit untersucht die Einstellorte von verschiedenartigen Teilchen sowie von 
kolloidalen und echten Lésungen im Feld stehender Ultraschallwellen. Spezifisch schwerere Teil- 
chen, die in einer Fliissigkeit schon vor Schalldurchgang suspendiert sind, werden in die 
Ebenen der Druckknoten getrieben. Bringt man hingegen eine wasserige Suspension in das schon 
bestehende Schallfeld, so stellen sich schwerere Teilchen (auch kolloidale) in die Druckbiuche, 
leichtere in die Druckknoten ein. Selbst echte Lésungen flieBen beim Einstrémen in die durch- 
schallte Fliissigkeit in die Bauche, wenn sie schwerer, und in die Knoten, wenn sie leichter sind 
als die durchschallte Fliissigkeit. Echte Lésungen gleicher Dichte wie die Grundfliissigkeit wan- 
dern sowohl in die Knoten, als auch in die Bauche. 


Summary 


This work concerns the location of particles in stationary ultrasonic waves. If heavy particles 
are suspended in a liquid before turning on the sound field they are driven to the pressure nodes. 
If, however, the field is already established and the particles then added, heavy ones are driven 
to the antinodes, light ones to the nodes. True solutions added to the irradiated liquid move 
toward the nodes if they are heavier, to the nodes if lighter than the irradiated liquid. If both 
have the same density, the solution may move either way. 


Sommaire 


Dans le présent travail, on détermine les régions ou se rassemblent dans un champ d’ultrasons 
stationnaires les particules de différente nature, ainsi que les solutions colloidales et les solutions 
véritables. Les particules spécifiquement les plus lourdes qui se trouvent déja en suspension dans 
un liquide avant le passage des ultrasons, se rassemblent dans les plans des neeuds de 
pression. Par contre, si on place dans un champ sonore déja établi une suspension aqueuse, les 
particules les plus lourdes (méme colloidales) se rassemblent aux ventres de pression, tandis que 
les particules les plus légéres se rassemblent aux nceuds de pression. Les solutions véritables 
elles-mémes, versées dans un liquide irradié aux ultrasons, s’écoulent vers les ventres si elles sont 
plus denses que le liquide irradié, et vers les nceuds si elles sont moins denses que ce liquide. Les 
solutions véritables qui ont méme densité que le liquide irradié s’écoulent aussi bien vers les 


neeuds que vers les ventres. 


Teilchen, die in ein Feld stehender Ultraschall- 
wellen in Fliissigkeiten gebracht werden oder 
schon vor Bestehen des.Feldes in der Fliissigkeit 
suspendiert sind, stellen sich in diesen in Ebenen 
im Abstand einer halben Schallwellenlange ein 
und vermitteln so ein Bild der stehenden Ultra- 
schallwelle in der Fliissigkeit. Diese Erscheinung 
wurde, von verschiedenen Gesichtspunkten aus- 
gehend, schon in zahlreichen Arbeiten beschrie- 
ben. Es wurden Teilchen, wie zum Beispiel Koks- 
staub, Quarzpulver oder rote Blutkérperchen, in 
Wasser gebracht, um durch ihr Sinken in den 
Knotenstellen der Schallwelle ein Bild derselben 
zu vermitteln (BoyLE, LEHMANN, SORENSEN, 
Doenon, BriaNncAnt u. a.) [1]. Gasblasen wurden 
in der gleichen Art verwendet (BoyLE und 
FROMANN, S6GLLNER) [2]. Auch bei Koagulations- 
erscheinungen wurde dhnliches beobachtet. So 
gibt Marinesco [3] an, daB bei Erzeugung einer 
stehenden Ultraschallwelle in einer konzentrier- 


ten waSrigen Lésung von Kongorot die gréberen 
Farbstoffteilchen in den Schwingungsbauchen 
koagulieren. SOLLNER und Bonpy [4] zeigen, da 
bei einer Fliissigkeit, die Teilchen enthalt, die 
teils schwerer und teils leichter als die Fliissigkeit 
sind, die leichteren in die Bewegungsknoten der 
stehenden Welle wandern und die schwereren in 
die Schwingungsbauche. Bei kolloidalen Zer- 
teilungen soll dies aber nicht mehr zutreffen. 

Ein Vergleich dieser zahlreichen Arbeiten zeigt 
wohl bei allen eine Einstellung der Teilchen im 
Abstand von einer halben Wellenlange, iiber den 
Einstellungsort der Teilchen (Knoten oder 
Bauche) herrscht jedoch keine Ubereinstim- 
mung. Dieser ist, wie die im folgenden beschrie- 
benen Versuche zeigen, durch die Art des Ein- 
bringens der Teilchen bedingt. 

Zur Untersuchung dieser Frage wurden die 
Teilchen im Schlierenbild der Knoten und Bauche 


beobachtet. Die Versuchsanordnung war folgende: 
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In einem mit Wasser gefiillten Glastrog wurde 
der Schallkopf eines normalen Ultraschallthera- 
piegerates mit vertikaler oder etwas geneigter 
Strahlflache angebracht (Frequenz 1000 kHz, 
Schallintensitat 2---3 W /cm?). Als Reflektor wurde 
eine Glasplatte im Abstand von einigen Zenti- 
metern dem Schallkopf gegeniiber gestellt. Die 
zur Feststellung des Einstellortes der Teilchen 
verwendete Schlierenanordnung ist in Abb. 1 
dargestellt. Die Beobachtung des Schlierenbildes 
erfolgte nach Ausblendung des direkten Strahlen- 
biindels durch Blende B auf einem Schirm S, auf 
dem gleichzeitig die in das Schallfeld gebrachten 
Teilchen scharf abgebildet wurden. 


Abb. 1. Schlierenanordnung zur Feststellung des Einstell- 
ortes der Teilchen. 


Die Versuche zur Einstellung der Teilchen 
wurden nun auf zwei verschiedene Arten ausge- 
fiihrt: 

1. Herstellung einer Suspension von Teilchen in 
Wasser mit nachtréglicher Erzeugung von 
stehenden Wellen in der Suspension. 

2. Einbringen einer waBrigen Aufschwemmung 
der Teilchen, in das schon in Wasser be- 


stehende Ultraschallfeld. 


Versuch nach 1: 


Von dem im Wasser erzeugten stehenden Ultra- 
schallfeld wurde ein Schlierenbild auf dem 
Schirm entworfen und auf diesem wurden die 
Stellen der Druckknoten der stehenden Ultra- 
schallwelle gekennzeichnet. Auf Vermeidung 
einer Parallaxe wurde dabei besonders geachtet. 
Nach Entfernung der Blende B wurde, ohne 
Ultraschalldurchgang, in der Fliissigkeit zwi- 
schen Schallkopf und Reflektor eine gleich- 
maBige Suspension von Weizenstarketeilchen im 
Wasser hergestellt und erst jetzt wieder die 
stehende Ultraschallwelle in der Flissigkeit er- 
zeugt. Die Teilchen stellten sich in Ebenen paral- 
lel zu Schallkopf und Reflektor in den Druck- 
knoten der stehenden Schallwelle ein. Bei Ver- 
wendung leichterer Teilchen (Emulsion von Ben- 
zol in Wasser) stellten sich die Benzoltrépfchen 
jedoch in den Druckbéuchen der stehenden Welle 


ein. 


Versuch nach 2: 


Wieder wurden die Stellen der Druckknoten 
mit Hilfe des Schlierenbildes auf dem Schirm an- 
gezeichnet. Dann wurde jedoch mit Hilfe einer 
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Pipette eine Aufschwemmung von Weizenstarke 
in das schon bestehende Feld gebracht. Dabei 
stellten sich nun die Weizenstarketeilchen in die 
Druckbduche der Schallwelle ein. Um den Ver- 
such mit leichten Teilchen durchzufiihren, wurde 
wieder eine grobe Emulsion von Benzol in Wasser 
hergestellt und in das Feld gebracht. Die Benzol- 
trépfchen begaben sich in die Druckknoten. Abb.2 
zeigt sowohl Weizenstarketeilchen als auch Ben- 
zoltrépfchen im Feld der stehenden Ultraschall- 
welle, wobei es sich bei der Abbildung um eine 
Photographie des Schirmbildes handelt. Die 
langen Striche sind Weizenstirketeilchen, die 
dazwischen befindlichen kugelférmigen Teilchen 
Benzoltrépfchen. 


Abb. 2. Weizenstiirketeilchen und Benzoltrépfchen im Feld 
der stehenden Ultraschallwelle. 


Die Tabelle I faBt die Ergebnisse iibersichtlich 


zusammen: 


Tabelle I 


Teilchen sind vor Be- 
stehen des Feldes in 
der Fliissigkeit gleich- 
maBig verteilt 


Teilchen werden in das 
schon bestehende Feld 
gebracht 


leichtere Teilchen Druckbiauche Druckknoten 


Druckbiauche 


schwerere Teilchen Druckknoten 


Das Bestehen der beiden dargelegten Gleich- 
gewichtslagen lieB sich auch in folgender Weise 
zeigen: Stért man etwa durch vorsichtiges Durch- 
fahren des Feldes mit einem Glasstab die Gleich- 
gewichtslage 1, so verlassen die in den Knoten 
eingestellten Weizenstarketeilchen teilweise ihre 
Platze und wandern in die Druckbauche. 
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Einstellung von Liésungen 


Die folgenden Versuche zeigen, daB sich im 
Gegensatz zu den bisherigen Erfahrungen auch 
kolloidal, ja sogar echt geléste Teilchen, aller- 
dings zusammen mit ihrem Lésungsmittel, ein- 
stellen. Bei diesen Versuchen wurde in die bisher 
verwendete Anordnung mittels einer Pipette ein 
Tropfen einer Kongorotlésung in das schon in 
Wasser, bestehende Ultraschallfeld eingebracht. 
Dabei bilden sich sofort rote Streifen mit dem 
Abstand einer halben Wellenlange. Diese Streifen 
werden durch Diffusion langsam unscharf und 
verschwinden schlieBlich trotz des bestehenden 


Feldes. 


Wird ein Tropfen einer echten Lésung auf die 
gleiche Weise wie oben eingefihrt, so tritt die 
gleiche Streifenbildung auf. Als echte Lésungen 
wurden Kaliumbichromat, Kupferchlorid, Eisen- 
chlorid und Kaliumpermanganat verwendet; die- 
ses zeigt wegen seiner starken Farbung den 
Effekt besonders schén. 

Da es sich bei den echt gelésten Teilchen 
keinesfalls um die Einstellung der Teilchen selbst, 
sondern nur um die der Lésung als Ganzes han- 
deln kann, wurde nach dem fiir die Einstellung 
der Lisung maBgebenden Merkmal gesucht und 
schlieBlich gefunden, dafi das Vorhandensein 
eines Dichteunterschiedes zwischen eingebrachter 
Lésung und dem Wasser wesentlich ist. 

Um den Einflu8 des Dichteunterschiedes nach- 
zuweisen, wurden wieder die Einstellungsorte 
der eingebrachten Lésung im Schlierenbild in der 
bekannten Versuchsanordnung beobachtet. Es 
wurde zunichst die Einstellung von Kalium- 
permanganatlésungen verschiedener Konzentra- 
tion in Zuckerlésungen verschiedenen spezi- 
fischen Gewichtes gepriift. Die Lésungen wur- 
den wieder in das schon bestehende Ultraschall- 
feld mittels einer Pipette eingebracht, mit dem 
Ergebnis, da jene Lésung, die leichter ist als 
die Zuckerlésung, in den Druckknoten aufwarts 
strémt, wahrend die schwerere Lésung in den 
Druckbéuchen sinkt. Von besonderem Interesse 
war nun das Einbringen einer Kaliumperman- 
ganatlésung von gleicher Dichte wie die Zucker- 
lésung. Dies fiihrte zur gleichzeitigen Einstellung 
der Lésung sowohl in den Druckbauchen als auchin 
den Druckknoten, somit im Abstand einer Viertel- 
Wellenlange, wie die verschiedenen Streifen in 
Abb. 3. 

Durch kleine Variationen der Dichte treten 
einmal die einen Streifen und einmal die anderen 
starker hervor. Quantitative Untersuchungen er- 
gaben, daS beim Abweichen der Dichte der ein- 
gebrachten Liésung von der der Grundlésung um 
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Abb. 3. ,,4/4-Streifen** bei gleicher Dichte der Fliissigkei- 


ten. 


1%, die beiden Streifensysteme gerade noch er- 
kennbar sind, bei gré®Beren Abweichungen aber 
nur eines. Das Auftreten dieser ,,j/4-Streifen‘ 
ist also ein sehr empfindliches Kennzeichen fiir 
die Dichtegleichheit der beiden Lésungen. Wird 
das Feld abgeschaltet, solange noch (/4-Streifen 
vorhanden sind, so sieht man die in den Druck- 
bauchen gebildeten Streifen entsprechend ihres 
Dichteunterschiedes gegeniiber der Zuckerlisung 
sinken, die anderen hingegen aufsteigen. 


Zusammenfassend ergeben sich fiir Lésungen 
die Einstellorte nach Tabelle II. 


Tabelle II 
Eingebrachte kolloidale oder echte Lésung a 
gegeniiber der vom Ultraschall durchsetzten Einstellort 
Fliissigkeit : 
leichter Druckknoten 
schwerer Druckbiauche 


Druckknoten und 
Druckbauche 


gleiche Dichte 


AbschlieBend sei noch ein Versuch angegeben, 
der die Einstellung von Lésungen in besonders 
einfacher und klarer Weise zeigt. Zur Beobach- 
tung dient eine kleine Kammer, die auf einfache 
Weise hergestellt werden kann. Der Raum iiber 
einer Glasplatte wurde mit Glasplattchen auf drei 
Seiten so abgeschlossen, da eine Kammer der 
GriéBe 4 x 4 cm? und der Héhe 1 mm entstand, 
die auf der vierten Seite durch eine diinne 
Gummimembran abgeschlossen wurde. Von hier 
konnte der Ultraschall zugeleitet werden. Die 
Zufiithrung des Ultraschalles erfolgte in der Weise, 
daB an das Tischchen, auf dem die Kammer lag, 
ein GefaB mit Wasser gestellt wurde, in das der 
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Abb. 4. Einbringen von Kaliumper- 
manganatlésung (dunkle 
Streifen) in Wasser. 


Schallkopf zwecks Kihlung zur Halfte ein- 
tauchte; die Kammer wurde mit der Membran 
an den Schallkopf angeschoben und zur besseren 
Schalliibertragung wurde zwischen den Schall- 
kopf und die Gummimembran ein Wassertropfen 
gebracht. (Diese Kammer eignet sich auch fir 
Ultraschallversuche mit gleichzeitiger Beobach- 
tung im Mikroskop.) 

Fullt man die Kammer mit Wasser, erzeugt 
dann eine stehende Ultraschallwelle und gibt 
einen Tropfen Kaliumpermanganatlésung darauf, 
so bilden sich rote Streifen (Abb. 4); erzeugt man 
jedoch das Schallfeld in einer Kaliumpermanga- 
natlésung und gibt einen Wassertropfen darauf, 
so sieht man helle Streifen auf rotem Grunde 
(Abb. 5). Beide Streifensysteme nebeneinander 
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Abb. 6. Streifenbildung durch Ein- 
bringen von Wasser (links) 
und konzentrierter Kalium- 
permanganatlésung (rechts) 
in Kaliumpermanganatlé- 
sung maBiger Konzentration. 


bilden sich, wenn man in eine mabig konzentrierte 
Kaliumpermanganatlésung von der einen Seite 
eine konzentrierte Lésung und von der anderen 
Seite Wasser tropft. Hierbei entstehen, wie Abb. 6 
zeigt, auf der einen Seite (rechts in der Abbildung) 
dunkle Streifen auf maBig dunklem Grunde und 
auf der anderen Seite helle Streifen auf demselben - 
Grunde. Die Streifen sind dabei um eine Viertel- 
Wellenlange gegeneinander verschoben. 

(Eingegangen am 2. September 1954.) 
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ETUDE STATISTIQUE 
DE LA SENSIBILITE AUDITIVE DIFFERENTIELLE D’INTENSITE 


par R. CHOCHOLLE 


Groupe de Recherches de Physiologie des Sensations, Laboratoire de Neurophysiologie du Collége de France 


Sommaire 


Les courbes de réponse 4 la variation d’intensité ont été étudiées a partir de divers niveaux 
et a trois fréquences. Ces courbes, de forme sigmoide, ont une pente raide; elles sont d’autant plus 
redressées que l’intensité est plus élevée; aux fortes intensités, elles deviennent asymétriques. Les 
seuils différentiels semblent indépendants de la fréquence; par contre, ils varient avec l’intensité, 


comme on le sait déja. 


Zusammenfassung 


Die Hérempfindlichkeit fiir Intensititsunterschiede wurde fiir verschiedene Bezugspegel bei 
drei Frequenzen untersucht. Die Kurven haben die Form einer Sigmoide und verlaufen sehr steil, 
und zwar um so steiler, je héher die Intensitit ist. Bei sehr groBen Intensititen werden sie unsym- 
metrisch. Die Unterschiedsschwellen scheinen unabhangig von der Frequenz zu sein, andern sich 


aber, wie schon bekannt war, mit der Intensitit. 


Summary 


Response curves of the ear to intensity variation were studied, relatively to different intensity 
levels, at three frequencies. These curves, of sigmoid form, exhibit a steep slope; the louder the 
tone, the steeper the slope; at high intensities, they are asymmetrical. Differential thresholds seem 
not to depend on frequency, but on the other hand, they vary with intensity. 
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1. Introduction 


Les essais faits antérieurement sur les seuils 
auditifs absolus m’avaient montré, entre autres 
résultats [1], [2], que les fluctuations des valeurs 
liminaires sont faibles, et, d’autre part, que les 
courbes cumulatives de la fréquence des réponses 
affirmatives en fonction de l’intensité du stimulus 
(plus couramment appelées courbes de réponse) 
ont une allure sigmoide fortement redressée. 

J’ai repris une étude statistique analogue sur 
les seuils auditifs différentiels d’intensité; j’ai 
étudié, par comparaison, la répartition du pour- 
centage de réponses affirmatives en fonction de 
l’amplitude de la variation d’intensité, et ceci a 
différents niveaux et différentes fréquences. 


2. Technique des mesures 
a) Appareillage 
La Fig. 1 donne un schéma-du montage utilisé. 
Les lignes étaient asymétriques et les impédances 


correctement adaptées, dans un sens comme dans 
VYautre; Vécouteur était un écouteur électro- 


dynamique. 
Attenuateur 
1 


FRE Mee 


Chambre insonore 


Oscillateur [| 


Trembleur 


Fig. 1. Schéma de l’appareillage utilisé. 


Une clé-inverseuse double a bascule (B), fixée 
prés du sujet dans la chambre insonore, per- 
mettait, simultanément, de mettre pour un temps 
hors du circuit l’atténuateur 2, et de fermer le 
circuit en dehors de cet atténuateur; les déplace- 
ments de quelques gouttes de mercure dans de 
petits godets a vide permettaient ce passage d’un 
circuit a l’autre. 

Normalement, l’inverseur se trouvait dans la 
position a; les atténuateurs | et 2 étaient, done, 
dans le cireuit; quand on appuyait sur le bouton 
de Vinverseur B (position b), l’atténuateur 2 
était mis hors du circuit, et il en résultait, évidem- 
ment, une augmentation temporaire de l’inten- 
sité; dés qu’on relachait le doigt, l’inverseur 
retournait de lui-méme A sa position initiale, et 
Vintensité revenait ainsi a sa valeur primitive. 
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Cette clé-inverseuse permettait, donc, de pro- 
duire des accroissements d’intensité de durée 
moyenne. D’autre part, au cours de l’inversion 
des circuits, la fermeture de l’un avait lieu trés 
légérement avant l’ouverture de l’autre; ainsi, 
malgré la brusquerie de la variation, il y avait un 
palier intermédiaire, d’oti une réduction des 
transitoires. D’une fagon générale, d’ailleurs, le 
dispositif utilisé a l’avantage de réduire consi- 
dérablement les transitoires (contacts de mercure 
dans le vide, butées en caoutchouc mousse, etc.) ; 
c’est en raison de ce fait que ce dispositif a été 
utilisé de préférence 4 d’autres, comme de nom- 
breux et longs essais l’ont prouvé. 

Une clé ordinaire (clé A) permettait de ne 
fermer le circuit général sur l’écouteur que le 
temps de chaque essai, ceci afin de ne pas fatiguer 
le sujet. 

L’opérateur tenait constante la somme des 
atténuations en décibels apportées par les deux 
atténuateurs dans la position a de la clé-inver- 
seuse, quelle que soit l’atténuation apportée par 
V’atténuateur Z; ainsi, les différents accroisse- 
ments temporaires d’intensité étaient toujours 
obtenus relativement a un niveau initia] toujours 
identique. Si le niveau initial choisi était N (dB) 
et les accroissements d’intensité désirés AN, AN’, 
AN’, --- (dB), Vatténuateur 2 était fixé a AN, 
AN’, AN", --- (dB), et Patténuateur 1 4 N—AN, 
N—AN', N—AN’,.:-- (dB). Le stimulus, aprés 
Vinversion du circuit, devenait, donc, N+ AN, 
N+AN’, N+AN", --- (dB). 

Ce montage a l’avantage d’éviter tout circuit 
secondaire dans les réseaux d’affaiblissement des 
atténuateurs (ce qui se produirait si on ne mettait 
pas entiérement hors du circuit l’atténuateur 2); 
et on peut étre certain de la précision des varia- 
tions introduites (a l’inverse d’un montage par 
substitution). 

Le sujet avait 4 sa disposition une clé de 
réponse insérée sur le circuit d’un trembleur 
situé a l’extérieur de la chambre insonore. 


b) Méthode des mesures 


La méthode des stimuli constants a été em- 
ployée dans ces mesures; mais, comme dans 
toutes les expériences antérieures, elle avait été 
modifiée de fagon que le sujet collabore active- 
ment aux essais [1], [2]. 

Un signal lumineux discret prévenait le sujet 
d’avoir a faire trois essais successifs. Dés qu'il se 
sentait prés, il appuyait sur le bouton de la clé 
A, et, aprés un temps d’au moins une seconde, il 
appuyait sur le bouton de la clé-inverseuse B, 
tout en maintenant le contact sur le bouton de la 
clé A (cette manceuvre était facilitée par la posi- 
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tion des deux clés: étant fixées céte a céte, on 
pouvait appuyer simultanément sur l’une avec 
l’index, et sur l’autre avec le pouce ou le majeur 
de la méme main); au bout d’une seconde au 
moins, le sujet relachait les deux clés. Ainsi, aprés 
étre parti du silence, il avait entendu d’abord 
le son au niveau initial N pendant un temps 
suffisant, puis, sans aucun intervalle, le méme 
son a un niveau N+ AN pendant un temps suffi- 
sant, lui aussi; enfin, il retournait au silence. Ce 
premier essai était suivi d’un repos de quelques 
secondes, puis d’un nouvel essai, et, avec le 
méme processus, d’un troisiéme. 

Quand le sujet avait terminé les trois essais 
demandés, il indiquait au moyen d’un signal 
convenu le nombre d’essais au cours desquels il 
avait percu une variation d’intensité. 

L’opérateur modifiait ensuite la valeur de 
Vaccroissement AN, tout en maintenant constant 
le niveau initial NV comme il a été dit antérieure- 
ment; le sujet faisait trois nouveaux essais, et 
ainsi de suite. 

I] était recommandé au sujet de laisser durer le 
son suffisamment longtemps, aussi bien avant 
qu’aprés avoir appuyé sur la clé-inverseuse B; 
on lui demandait aussi de ne s’attacher qu’a 
Vaspect intensité du stimulus avant et aprés 
entrée en action de la clé B, et ceci d’une 
fagon globale, sous forme d’une véritable com- 
paraison, en évitant soigneusement de préter 
attention a ce qu’il pourrait ressentir au moment 
ou il appuie sur la clé B. 

Les essais étaient faits en suites croissantes et 
décroissantes alternées; cette méthode permet de 
compenser les sous-estimations et les surestima- 
tions provenant du sens des essais. Les séries 
étaient elles-mémes fractionnées, de fagon A ne 
pas étre trop longues: on commengait par une 
suite décroissante de variations AN d’intensité a 
partir d’un méme niveau initial N, suivie d’une 
suite croissante aux mémes valeurs; les bonds 
étaient de 0,3 en 0,3 dB; cette demi-série était 
suivie immédiatement d’une autre demi-série 
identique a la précédente, avec la seule différence 
que tous les accroissements AN testés étaient 
supérieurs de 0,15 dB aux précédents. Ainsi, au 
total, par rapport 4 un méme niveau initial N, 
les divers accroissements testés différaient entre 
eux de 0,15 dB. 

Tous les trois ou quatre accroissements d’in- 
tensité testés, on présentait au sujet un accroisse- 
ment de 10 dB plus grand, a la méme fréquence 
et a partir du méme niveau initial N. 

A chaque niveau initial utilisé, les accroisse- 
ments AN testés étaient répartis de part et 
d’autre de la valeur du seuil différentiel; 12 a 14 
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valeurs furent testées aux forts niveaux d’inten- 
sité, et 16 A 18 valeurs aux faibles niveaux. A 
chaque valeur testée, 100 essais furent effectués. 
Chaque série d’essais durait 15 4 30 minutes; elle 
était suivie d’un repos de 15 a 20 minutes. 

Avant chaque série, on procédait 4 une déter- 
mination rapide du seuil absolu a la fréquence 
utilisée, par la méthode des stimuli constants. 
Cette détermination était rendue nécessaire pour 
établir le niveau initial: la valeur de ce niveau, N, 
était, en effet, fixée par rapport au seuil absolu; 
de cette facon, les variations du seuil absolu 
d’une série a |’autre n’intervenaient pas dans les 
résultats. Les niveaux initiaux JV utilisés furent: 
1, 3, 6, 10 dB au-dessus du seuil, et, au-dela de 
cette valeur, tous les niveaux possibles de 10 en 
10 dB. Les sujets oni travaillé 4 trois fréquences: 
200, 1000 et 10000 Hz. 

J’ai fait appel a quatre sujets bien entrainés 
et A un 5¢ sujet peu entrainé. Les quatre premiers 
sujets ont travaillé plusieurs heures par jour, 
pendant de longs mois. 


3. Résultats 


Les figures 2 4 7 résument les résultats pour 
deux des sujets, a trois fréquences. 

Les abscisses sont les accroissements d’inten- 
sité AN en décibels, 4 partir du niveau initial N 
indiqué pour chaque courbe, et les ordonnées, 
les pourcentages de réponses affirmatives. 
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Fig. 2. Pourcentage de réponses affirmatives en fonction 
des accroissements d’intensité AN, évalués en déci- 
bels; chaque courbe correspond 4a un niveau NV 
d’intensité initialement a la variation; les niveaux 
auxquels chaque courbe a été obtenue sont indiqués 
a droite. Sujet I 4 200 Hz, en écoute monaurale. 


Fig. 3. Mémes explications que pour la Fig. 2; 
Sujet 1 a 1000 Hz. 
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Fig. 4. Mémes explications que pour la Fig. 2; 
Sujet I a 10000 Hz. 
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Fig. 5. Mémes explications que pour la Fig. 2; 
Sujet II a 200 Hz. 


Fig. 6. Mémes explications que pour la Fig. 2; 
Sujet II a 1000 Hz. 
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Fig. 7. Mémes explications que pour la Fig. 2; 
Sujet IT a 10 000 Hz. 


La forme des courbes de réponse a été déduite 
de celles obtenues dans une étude antérieure sur 
les seuils absolus monauraux et binauraux [1], [2]. 

En effet, la méme forme de courbes s’applique 
exactement aux essais faits aux niveaux initiaux 
N juxtaliminaires; aux niveaux initiaux N plus 
intenses, il suffit, semble-t-il, de modifier la 
pente, tout en conservant l’allure générale, pour 
avoir des courbes convenables; il suffit de choisir 
judicieusement le rapport de transformation. 

Aux intensités élevées cependant, les courbes 
de réponse sont nettement dissymétriques. 
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4. Discussion 


a) Forme des courbes de réponse 


Ainsi qu’il vient d’étre dit, la forme représenta- 
tive des courbes de réponse des seuils différen- 
tiels en dB aux niveaux initiaux N juxtali- 
minaires est en tous points la méme que celle 
des seuils absolus; ces deux séries de courbes sont 
superposables, pour un méme sujet. Ceci est 
d’autant plus évident dans mes expériences que 
des sujets furent communs aux séries d’essais 
sur les seuils absolus et les seuils différentiels. 

Aux intensités plus élevées, la forme générale 
subsiste; seule, la pente change. 

Ainsi, il n’y a pas de discontinuité entre sensi- 
bilité absolue et sensibilité différentielle. On peut, 
donc, dire que la sensibilité absolue n’est qu’un 
cas particulier de la sensibilité différentielle, dans 
lequel le stimulus au niveau initial N est le 
silence, N étant, alors, au plus égal a 0; on 
retrouve ainsi ce que j’avais déja montré dans 
une étude antérieure sur les temps de réaction [3]. 

Dans aucun cas, les points ne peuvent étre 
joints par une droite, comme l’avaient fait divers 
auteurs: STEVENS, Morcan et VOLKMANN [4], 
[5], von BéKésy [6], a partir de trop peu de 
points, il est vrai (trois ou quatre points, le plus 
souvent). 

La forme sigmoide des courbes de réponse s’op- 
pose formellement a l’hypothése émise par ces 
auteurs; d’aprés eux, il existerait une échelle de 
quanta nerveux, de neuroquanta; le déclenche- 
ment de N—1, N, N+1, ---, neuroquanta, 4 un 
certain niveau du stimulus, dépendrait du niveau 
atteint par un processus qui varierait plus ou 
moins dans le temps; mais si les chances d’appari- 
tion sont égales pour tous les neuroquanta N—1, 
N, N+1,.---, la courbe de réponse devrait étre 
rectilinéaire. Or, ceci ne se vérifie pas. D’autre 
part, il devrait résulter de cela que, pour passer 
d’une extrémité de la courbe a l’autre extrémité, 
on devrait étre obligé de doubler la variation 
d’intensité; or, ceci ne se vérifie pas non plus. 

Les résultats précédents sont en accord, par 
contre, avec ceux d’autres auteurs [7], [8]. 


b) Forme asymétrique des courbes de réponse 


Aux niveaux élevés d’intensité, la forme des 
courbes de réponse devient nettement asymétri- 
que, et cela d’autant plus que le niveau est plus 
élevé. Cette asymétrie est encore plus nette si on 
traduit les accroissements d’intensité en valeurs 
absolues. Ceci est en accord, d’ailleurs, avec les 
résultats de Monrcomery [7]. 

Ce fait ne doit pas nous étonner. En effet, le 
seuil différentiel est trés faible A ces niveaux, si 
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bien que l’évolution de la portion des courbes de 
réponse au-dessous des 50% de réponses affirma- 
tives ne peut plus se dérouler normalement, tan- 
dis qu’elle peut se dérouler normalement dans la 
moitié supérieure. 


c) Pente des courbes de réponse 


La pente est trés redressée; elle est, d’ailleurs, 
d’autant plus redressée que le niveau NV auquel 
on fait les essais est plus élevé, comme [avait 
déja constaté MontcomeEry [7]. 

Aux niveaux juxtaliminaires d’intensité (1dB 
au-dessus du seuil) la pente est, comme on I’a vu, 
sensiblement identique 4 celle des courbes de 
réponse pour la sensibilité 4 Vapparition d’un 
stimulus sonore, pour un méme sujet et a la 
méme fréquence; il suffit de faire varier les 
accroissements de l’intensité AN de 1,2 4 1,5 dB 
pour passer de 0 a 100 % de réponses affirmatives. 

A des niveaux N de 10 4 20 dB au-dessus du 
seuil la pente est, déja, bien plus raide: on passe, 
alors, de 0 a 100% de réponses affirmatives en 
variant les accroissements de l’intensité de 0,5 a 
0,7 dB. 

Aux niveaux supérieurs 4 40 dB au-dessus du 
seuil, la pente est encore plus raide; aux niveaux 
élevés, toutes les courbes de réponse évoluent dans 
une marge de 0,3 a 0,4 dB. 

Ce redressement progressif avec le niveau N de 
Vintensité permet de répondre 4 l’objection sui- 
vante: les transitoires, bien qu’inaudibles, n’in- 
terviendraient-ils pas dans une certaine mesure ? 
Il est bien certain que, malgré les précautions 
prises, on n’est pas a l’abri de transitoires; il 
peut encore s’en créer au niveau de |’écouteur ou 
de loreille; bien qu’inaudibles, et malgré les 
consignes données aux sujets, on pourrait se 
demander s’ils ne pourraient pas, cependant, jouer 
un réle. Mais, s’il en était ainsi, toutes les courbes 
de réponse a laccroissement d’intensité ne 
seraient, alors, que des courbes de réponse A 
Papparition des transitoires, c’est-a-dire a l’appa- 
rition d’un stimulus sonore; et toutes les courbes 
de réponse devraient avoir la méme pente, celle 
des courbes obtenues pour les seuils absolus [1], 
[2], celle encore des courbes de réponse obtenues 
aux niveaux juxtaliminaires. Or, il n’en est rien, 
méme aux niveaux N d’intensité ot, pourtant, 
les chances d’avoir des transitoires sont beaucoup 
plus grandes. 

Tout au plus, a ces niveaux élevés, peut-on se 
demander si la disposition, aberrante A premiére 
vue, de quelques points en dehors des courbes 
théoriques, n’est pas due a une certaine inter- 
vention de transitoires. A ces niveaux élevés, en 
effet, dans la partie supérieure des courbes de 
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réponse, on voit des points s’éloigner progressive- 
ment des courbes théoriques; or, plus la varia- 
tion introduite est grande, plus importants 
doivent étre les transitoires. Mais, méme s’il en 
était ainsi, l’intervention des transitoires dans les 
résultats serait trés limitée, puisqu’elle ne de- 
viendrait sensible qu’ l’extrémité supérieure des 
courbes de réponse. 

Quelle est la raison de cette forme sigmoide 
des courbes de réponse ? Ces courbes rappellent 
les courbes de répartition statistique: répartition 
de Gauss, répartition de Poisson reposant sur un 
grand nombre d’événements. Existe-t-il des élé- 
ments qui permettraient une telle interprétation ? 
Evidemment, il pourrait y avoir des fluctuations 
dans le temps a différents niveaux, que ce soit 
dans le stimulus, au niveau des éléments senso- 
riels, dans les voies nerveuses ou dans les centres 
[1], [2]. Il est impossible de répondre positive- 
ment ou négativement a la question posée. Ce- 
pendant, on ne voit pas trés bien comment on 
pourrait expliquer, dans ce cas, pourquoi il est 
impossible de tracer une courbe théorique de 
répartition statistique qui puisse s’accorder 
exactement avec les points expérimentaux; 
méme dans les meilleurs des cas, si on arrive a 
une coincidence excellente sur une bonne partie 
du parcours, des divergences nettes, souvent 
importantes, apparaissent toujours aux extré- 
mités des courbes. 


d) Variation des seuils différentiels avec l’intensité 


J’ai pris comme valeurs des seuils différentiels 
(SD) d’intensité aux divers niveaux N les 
accroissements d’intensité AN qui permettent 
d’avoir 50% de réponses affirmatives. Les figures 
8 a 10 résument les résultats; ces courbes mon- 
trent les variations des seuils différentiels SD, 
évalués en décibels, en fonction du niveau WN de 
Vintensité. On retrouve ce qu’on savait déja, a 
savoir la décroissance progressive des seuils 
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Fig. 8. Variation des seuils différentiels SD, évalués en dé- 
cibels, en fonction du niveau initial de intensité N, 
indiqué en décibels au-dessus du seuil. Résultats 
individuels pour plusieurs sujets, en écoute monau- 
rale, 4 200 Hz. 
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Fig. 10. Mémes explications que pour la Fig. 8 a 10000 Hz. 


différentiels, au fur et 4 mesure que l’intensité 
croit. 

Il semble méme que cette décroissance con- 
tinue encore, bien qu’extrémement faiblement, 
aux niveaux élevés. Or, comme dans tous les 
phénoménes biologiques, on devrait atteindre un 
plateau. Ce plateau n’a encore jamais été trouvé 
par aucun chercheur dans |’audition. Riesz lui- 
méme [9], qui avait travaillé 4 1000 Hz jusqu’a 
90 dB, n’a pas trouvé de plateau; tout au plus, 
les seuils trouvés par lui semblent-ils tendre vers 
une valeur limite. 

Doit-on conclure que les phénoménes auditifs 
se comporteraient autrement que tous les autres 
phénoménes biologiques, ce qui serait inconce- 
vable ? 

Evidemment, aux niveaux élevés, les points 
sont trés rapprochés et, parfois méme, imbriqués 
les uns dans les autres (voir figures 2 4 7); cepen- 
dant, dans l’ensemble, ils semblent bien s’aligner 
sur des courbes différentes, aux différents niveaux 
N choisis, et ces courbes paraissent bien s’éche- 
lonner; ceciapparait encore plus clairement quand 
on utilise de grandes échelles. Le tracé des cour- 
bes sur les figures tient compte de cette tendance. 

Cependant, la réduction du seuil différentiel, 
quand croit le niveau N de l’intensité, finit par 
étre extrémement minime, aux niveaux élevés: 
0,03 a 0,04 dB, quand on passe d’un niveau N 
4 un autre niveau N’ de 10 dB plus intense. Or, 


pour 100 essais, d’aprés la relation op = Vpqin 
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(p représentant la fréquence des réponses affirma- 
tives, q celle des réponses négatives, et n le 
nombre de mesures), l’erreur-type op devrait étre, 
en théorie, de 0,5% pour p= 0,5; d’aprés la 
pente des courbes auprés des 50% de réponses 
affirmatives, cette erreur-type représenterait 


+ 0,01 4 0,02 dB, environ. 


Peut-on, dans ces conditions, affirmer qu’on 
n’atteint jamais un plateau dans l’audition ? On 
ne le peut vraiment pas, d’autant plus que, 
comme on l’a vu, peut-étre pourrait-il y avoir 
une intervention de transitoires, bien que cer- 
tainement fort minime; d’autre part, a ces 
niveaux élevés, ne peut-il pas intervenir une 
autre modalité sensorielle, vibratoire par exem- 
ple? est-on encore bien sar des valeurs du stimu- 
lus avec une telle précision, a de telles intensités, 
et de sa qualité ? Enfin, l’effort de Vattention 
pour essayer de percevoir une variation est consi- 
dérable a ces niveaux élevés; ceci pourrait bien 
se traduire dans les résultats. 


Et d’ailleurs, ce plateau pourrait bien se trou- 
ver au-dela. Mon appareillage ne m’a pas permis, 
en effet, d’aller au-dela de 100 dB au-dessus du 
seuil a 1000 Hz; mais de 100 4 120 ou 130 dB, il 


y a une marge encore importante. 


Ce qu’on peut dire de certain, en tout cas, c’est 
que la variation du seuil différentiel est bien 
minime, a partir de 40 a 50 dB au-dessus du seuil: 
elle est, au plus, de 0,2 4 0,4dB, quand on passe 
de 40 (ou 50) dB 4 100 dB au-dessus du seuil. Aux 
niveaux élevés de l’intensité, le seuil différentiel 
est de 0,15 4 0,25 dB, suivant les sujets. Je re- 
trouve ainsi a peu prés les valeurs données par 
Riesz [9], mais obtenues par une autre méthode. 


A un niveau NV d’intensité de 20 dB au-dessus 
du seuil, le seuil différentiel est de 0,5 4 1,2 dB; 
ce n’est qu’en-dega seulement de 10dB au- 
dessus du seuil que la sensibilité différentielle 
devient vraiment moins bonne; aux niveaux 
juxtaliminaires, elle atteint 2,5 4 3 dB. 


D’une fagon générale, les seuils obtenus dans 
ce travail sont un peu plus faibles que ceux 
donnés par d’autres auteurs, par exemple RriEsz 
[9], Dimmick et Oxson [10], Mrirer [11], 
Liiscuer et Zwisiocki [12], Harris [13], Mont- 
GOMERY [7], Corso [8]. 

Les différences constatées entre les résultats 
présents et ceux donnés dans d’autres travaux 
s’expliquent, probablement, par des différences 
dans les méthodes utilisées; Dimmicx [13], 
Montcomery [7], Vont, d’ailleurs, montré, en 
comparant les résultats obtenus avec différentes 
méthodes, ou en variant les consignes données 


aux sujets. 
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Les causes principales de ces différences doi- 
vent étre certainement recherchées dans lutilisa- 
tion de sujets entrainés, et de statistiques de 
résultats sur des individus, et non sur toute une 
population, comme cela a été le cas dans beau- 
coup de ces travaux. D’autre part, dans les 
expériences présentes, les sujets collaboraient 
activement aux essais, comme on l’a vu; ils 
étaient, de plus, dans le meilleur état possible au 
moment des essais; enfin, leurs réponses étaient 
d’autant plus précises quwils ne répondaient 
qu’aprés avoir fait trois essais successifs. 

Il y a, en plus, le plus souvent, des différences 
importantes entre la méthode de mesures utilisée 
ici et celles employées par d’autres auteurs. Cer- 
tains ont utilisé des stimuli brefs, mais sans se 
préoccuper des transitoires; d’autres se sont 
contentés de déterminer la plus petite variation 
audible, le sujet indiquant quand il entend 
apparaitre une variation, ou au contraire, quand 
celle-ci semble disparaitre. D’autres encore ont 
intercalé un intervalle de silence entre la cessa- 
tion du premier signal et l’apparition du deuxiéme 
(Upron et Hoxtway, [15]); or, la durée de cet 
intervalle intervient. 

Toutes ces derniéres années d’ailleurs, la plu- 
part des recherches entreprises sur ce sujet ont 
utilisé un autre type de stimulus; afin d’éviter 
les transitoires, les auteurs utilisérent des signaux 
modulés; parmi les travaux les plus connus, on 
peut citer, par exemple, Riesz [9], Liscuer et 
Zwis.ocki [12]. A un rythme optimum de modu- 
lation, on recherche la dénivellation qui permet 
d’avoir encore une sensation de fluctuation, ou 
la dénivellation 4 laquelle cette sensation com- 
mence a apparaitre. C’est, en plus, une méthode 
qui demande un minimum au sujet, elle est 
facile A utiliser; c’est, d’ailleurs, la méthode 
employée en audiométrie pour détecter la pré- 
sence éventuelle d’un recrutement bilatéral dans 
certains types de surdités (méthode dite du test 
de LiscHer). Mais est-on bien sir de mesurer, 
par cette méthode, véritablement un seuil diffé- 
rentiel ? On mesure bien un seuil de sensation 
papillotante; mais quelle est la relation entre le 
seuil différentiel et le seuil de sensation papillo- 
tante, a supposer méme qu'il y en ait une, 
comme on a trop tendance a l’admettre a pri- 
ori ? 

I] ne faut pas oublier, en effet, que ce dont nous 
prenons connaissance dans nos expériences, ce ne 
sont pas les effets au niveau de la cochlée; c’est 
au contraire, le plus souvent, le résultat de juge- 
ments de la part des sujets sur ce qu’ils ont pergu 
ou cru percevoir. Les sujets jugent d’aprés ce 
qu’a recu le cortex, aprés toutes les pertes, les 
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intégrations, les associations, qui ont pu se pro- 
duire A tous les niveaux de la voie sensorielle. 

Or, les jugements ne sont pas les mémes quand 
il faut juger de l’existence ou non d’un papillote- 
ment, ou de la présence ou non d’une variation; 
d’autre part, les indices sur lesquels repose le 
jugement auprés du seuil différentiel peuvent 
différer dans l’un et dans l’autre cas; on sait qu’un 
stimulus modulé entraine une rythmicité a4 une 
cadence identique au niveau de certains éléments 
corticaux; ceci ne pourrait-il pas devenir l’indice 
sur lequel s’appuie notre jugement, dans le cas 
d’un stimulus modulé ? 

Le seuil différentiel relatif, évalué en pression 
acoustique, passe de 40% environ aux niveaux 
juxtaliminaires 4 10 ou 15% a 20 dB, et mest 
plus que de 2 4 3% aux niveaux élevés; traduit 
en puissance acoustique, le seuil différentiel rela- 
tif passe, dans les mémes conditions, de 100 a 
5 %. Ces chiffres peuvent étre traduits en valeurs 
absolues. En valeurs absolues de pression acous- 
tique, le seuil varie relativement peu jusqu’a un 
niveau de 20 dB au-dessus du seuil: sa valeur 
absolue croit seulement de deux a trois fois; mais 
au-dela, le seuil, évalué en valeurs absolues de 
pression acoustique, croit rapidement; il arrive, 
a un niveau de 100 dB par exemple, a corres- 
pondre a une variation de pression acoustique 
2000 ou 3000 fois plus ample; traduit en valeurs 
absolues de pression acoustique, le seuil différen- 
tiel est des centaines de millions de fois plus 


ample A 100 dB qu’a 1 dB. 


e) Les échelons différentiels d’intensité 


On peut intégrer les valeurs des seuils différen- 
tiels, comme |’a fait Riesz [15]; on obtient, alors, 
le nombre n d’échelons différentiels discrimi- 
nables, les uns a la suite des autres, entre NV et 
N’ dB. J’ai procédé plus simplement, comme je 
Vavais déja fait dans un travail antérieur [3]; 
il suffit d’évaluer de proche en proche le nombre 
d’échelons dans des intervalles limités, d’aprés 
les valeurs en décibels des seuils différentiels dans 
chacun de ces intervalles; cette méthode donne 
des résultats plus proches de la réalité, comme je 
lai constaté, en essayant de résoudre le probléme 
inverse: connaissant les courbes reliant le nombre 
n d’échelons différentiels discriminables d’inten- 
sité, en déduire le seuil différentiel 4 un niveau IV 
quelconque. 

Le résultat final apparait dans les courbes des 
figures 11 a 13. 

Ces courbes sont pratiquement semblables a 
celles qu’avait données Riesz. Comme les seuils 
différentiels sont légérement plus faibles, il en 
résulte, évidemment, qu’a chaque niveau, le 
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200 


80dB 


Fig. 11. Nombre n d’échelons différentiels en fonction du 
niveau initial NV de l’intensité, NV étant évalué en 
dB au-dessus du seuil. Résultats individuels pour 
plusieurs sujets, en écoute monaurale 4 200 Hz. 


0 20 40 60 


80dB 


Fig. 13. Mémes explications que pour la Fig. 11 4 10000 Hz. 


nombre d’échelons différentiels est un peu plus 
grand sur les courbes présentées ici. 

On sait que la question n’est pas encore entiére- 
ment résolue de savoir si cette échelle est une 
véritable échelle de sensation, si, donc, elle se con- 
fond avec ce qui serait une «échelle de sonie», ou si 
ces deux fonctions sont différentes l’une del’autre. 


f) Indépendance des seuils différentiels d’intensité 
d’avec la fréquence 


Les seuils différentiels ne semblent pas varier 
avec la fréquence; 4 niveaux identiques d’inten- 
sité par rapport au seuil absolu, les seuils différen- 
tiels paraissent du méme ordre; il y aurait, ce- 
pendant, peut-étre une tendance vers des valeurs 
légérement plus fortes aux fréquences moyennes. 

Ceci est 4 rapprocher de l’absence de différences 
dans l’amélioration du seuil en écoute binaurale 
par rapport a l’écoute monaurale, quelle que soit 
la fréquence (CHOCHOLLE [1]). 

Pourtant, dans les expériences de Riesz [9], 
la variation avec la fréquence paraissait nette; 
les seuils différentiels paraissaient notablement 
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plus faibles aux fréquences moyennes, surtout 
aux faibles niveaux d’intensité. 

Mais, plus récemment, LiiscHer et ZWISLOCKI 
[12] n’ont pas trouvé, eux non plus, de différences 
avec la fréquence, du moins aux niveaux élevés 
d’intensité; aux niveaux moyens, il y aurait, 
peut-étre, une différence, mais elle serait bien 
minime, et, en tout cas, dans le sens indiqué par 
RIEsz. 

Ici encore, la technique employée doit étre 
responsable des différences. Moi-méme d’ailleurs, 
dans des mesures faites par une autre méthode 
et en rapport avec d’autres expériences [3], 
j avais retrouvé des résultats comparables 4 ceux 
de Rixsz. Il faut dire aussi que les moyens utilisés 
alors étaient plus rudimentaires et la technique 
moins élaborée. 


5. Conclusions 


Cette étude systématique montre que notre 
sensibilité a l’accroissement de lintensité est 
extrémement fine dans le domaine auditif; de 
plus, si les conditions sont excellentes, on trouve 
une trés grande précision dans les réponses, 
puisqu’aux niveaux élevés, il suffit de variations 
de 2 4 3% dans le niveau de pression acoustique 
du stimulus pour passer de 0 a 100 % de réponses 
affirmatives. 

Toutefois, la traduction des résultats en 
valeurs absolues montre qu’en réalité, ceci corres- 
pond a une variation fort importante en valeurs 
absolues (méme s’il s’agit d’évaluations en varia- 
tions de pression acoustique), aux niveaux élevés. 
Ceci pose une question sur la nature du stimulus 


au niveau de la cochlée. (Regu le 25 octobre 1954.) 
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MEASUREMENTS OF THERMAL DIFFUSION IN H,—He MIXTURES 


BY USING ULTRASONICS 


by A. VAN ITTERBEEK and J. NIHOUL 


Instituut voor Lage Temperaturen en Technische Physica, Leuven, Belgium 


Summary 

An ultrasonic interferometer has been used to measure the change of the concentration produced 
by thermal diffusion in H,-He mixtures. With this method the separation degree has been 
measured as a function of the concentration and the temperature. From these measurements 
we calculated the thermal diffusion factor x, thus we obtained «= 0:49. We measured the 
change of the concentration as a function of time in H,—N, mixtures by means of the velocity. 
We obtained a good agreement with CLARK JoNES and FurrRy’s theory (Rev. mod. Phys. 18 
[1946], 151). 


Sommaire 


Un interferométre acoustique a été utilisé pour mesurer la variation de la concentration pro- 
duite par la thermodiffusion dans des mélanges H.—He. Au moyen de cette méthode nous avons 
mesuré le degré de séparation en fonction de la concentration et de la température. A partir de 
ces mesures nous avons calculé le facteur de thermodiffusion «, ainsi nous avons obtenu « = 0,149. 
Nous avons mesuré la variation de la concentration en fonction du temps dans des mélanges 
H,—N, en mesurant la vitesse du son. Les résultats obtenus sont en accord avec la théorie de 
CLARK JONES et Furry (Rev. mod. Phys. 18 [1946], 151). 


Zusammenfassung 


Mit einem Ultraschallinterferometer wurde die durch Thermodiffusion hervorgerufene Kon- 
zentrationsinderung in H,—He-Gemischen gemessen. Auf diese Weise wurde das Entmischungs- 
verhiltnis als Funktion der Konzentration und der Temperatur bestimmt und daraus die ther- 
mische Diffusionskonstante « zu 0,149 berechnet. Durch Messung der Schallgeschwindigkeit be- 
stimmte man die zeitlichen Konzentrationsinderungen in H,—N,-Gemischen. Die gewonnenen 
Ergebnisse sind in guter Ubereinstimmung mit der Theorie von CLARK JoNES und Furry (Rev. 


mod. Phys. 18 [1946], 151). 


1. Introduction 


Two years ago, one of us together with ForREz 
and MAriEns [1] developed a method to measure 
the change of the concentration produced by 
thermal diffusion in gas mixtures with the use 
of an ultrasonic interferometer. At that time we 
applied the method to H,—N, mixtures. Since 
then the accuracy of the method has been 
increased and moreover the apparatus has 
undergone some tests. 

We have now applied the method to measure 
thermal diffusion in H,—-He mixtures as a 
function of concentration and temperature, the 
lowest temperature (T.) of the thermal diffusion 
apparatus being at the boiling point of liquid 
propane, -oxygen or -nitrogen respectively. The 
highest temperature T, corresponds to room 
temperature. As is well known, the change of the 
concentration is given by the equation: 


Ay = ky In(T,/T3) (1) 


with kp = xy,y2: y; and y, being the concentra- 
tions of the two components. In general it is 


found that the factor « is independent of tem- 
perature. In certain mixtures however a decrease 
of « at the lowest temperatures has been found. 
Three years ago one of us together with DE 
TROYER and RIETVELD [3] found a decrease of « 
down to liquid hydrogen temperatures. As our 
present method is more accurate than our earlier 
one we proposed to find out whether this decrease 
does not start already at liquid oxygen tempera- 
tures. From our measurements we calculated the 
thermodiffusion factor. The measurements could 
not be carried out down to liquid hydrogen 
temperatures for the reason that at those tem- 
peratures the velocity of sound becomes a fune- 
tion of pressure, but we hope, by adapting the 
method, to make it possible that the measure- 
ments can be made down to these low tempera- 
tures. 


2. Experimental apparatus and method 


The acoustical interferometer combined with 
the cryostat is represented in Fig. 1. It is on the 
same principle as that published in our earlier 
paper [1], to which we refer for the description. 
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To determine the change of the concentration, 
we make use of the following equation 


Gnidia mr a 
y+. Miyit Moye 


in which W represents the velocity of sound, cy 
and c, the specific heats, M, and M, the molar 
masses, R= 8.316 x 10’ erg, T absolute temper- 
ature. Our interferometer 
permits one to measure the 
length with an accuracy 
of +5 microns, so that 
the velocity of sound can 
be measured with a rela- 
tive error of 0.01%. But 
we hope to reach a still 
better accuracy in future. 

Another advantage of 
the method as compared 
with the earlier ones is that 
now it is also possible to 
measure the change of the 
concentration as a func- 
tion of time, besides the 
final concentration. Thus 
we were able to check the 
theory of CLARK JONES 
and Furry [4] by means 
of our earlier measure- 
ments on the H,—N, mix- 
tures as a function of pres- 
sure. 

The way of carrying out the measurements is 
as follows. The apparatus is first evacuated by 
means of a diffusion pump. Then the cryostat is 
filled with liquid oxygen or liquid nitrogen and 
the interferometer also is filled with hydrogen 
gas in order to precool the apparatus. After two 
hours the hydrogen is pumped away and the 
interferometer is filled with the mixture, which 
has been prepared before. The concentration is 
determined from the partial pressures and con- 
trolled by means of a velocity measurement in 
the mixture at room temperature. 

The thermodiffusion-velocity measurements 
are started immediately after the interferometer 
is filled with the mixture. The change of velocity 
is continuously measured as a function of time. 
From these measurements the variation of the 
concentration (Ay,) within the lowest part of the 
apparatus can be determined by means of Eq. (2). 
The total change of concentration Ay is cal- 
culated by means of Eqs. (3) and (4) 


1. Interferometer 
in cryostat. 


Fig. 


Ay = Ay, + Aye (3) 
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: Von ul: 
with Ay, =Ay, = V, : T, : (4) 


(Ay, is the change in the volume V, which is at 
room temperature 7’, and Ay, corresponds to the 
volume V, which is at the low temperature T,). 
For the ratio V,/V, we found experimentally 
18.5. The measurements were done at a frequency 


of 507.88 ke/s. 


3. Experimental results on H.-He mixtures 


Table I 


Concentra- 
tion 


Sa: TS a ed cm 


ea 
290.1 | 228 26.0 21.2 0.09 0.4 0.5 
290.1 | 228 26.0 33.8 0.06 0.8 0.9 
290.1 | 228 26.0 48.0 0.08 al 1.2 
290.1 | 228 26.0 53.4 0.07 1.0 1.1 
290.1 | 228 26.0 74.0 0.06 0.9 1.0 
291.2 | 90.1 | 27.0 11.0 0.2 1.4 1.6 
291.4} 90.1 | 28.0 32.6 0.6 3.5 4.1 
291.2} 90.1 | 28.4 34.7 0.6 3.5 4.1 
291.1} 90.1 | 26.0 54.6 0.7 3.7 4.7 
290.4] 90.1 | 27.0 67.2 0.6 3.3 3.9 
288.1) 90.1 | 28.0 82.2 0.3 2.0 2.3 
290.6 | 77.3 | 26.0 13.6 0.5 2.4 2.9 
290.2'|' 77.3 | 25.0 34.5 0.7 3.7 4.4 
290.6 | 77.3 | 26.0 52.3 0.9 4.3 Bed 
287.9 | 77.3 | 27.0 56.6 0.8 4.0 4.8 
289-1 || 77.3) || 27.0 58.6 0.8 4.0 4.8 
288.1 | 77.3 | 26.0 81.4 0.6 2.8 3.4 


The experimental results are given in TableI and 
graphically drawn in Fig. 2, 3 and 4 respectively. 


In Fig. 5 we have graphically drawn Ay as a 
function of In(T,/T;,) 
* | ‘T= 228K 
; A 
0 20 40 60 80  0%He 


——H) He—+ _ —+¥ 


Fig. 2. Ay as a function of y at 228°K. 
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Fig. 3. Ay as a function of y at 90.1°K. 
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Fig. 4. Ay as a function of y at 77.3°K. 
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Fig. 5. Ay as a function of In(T,/T,) for y = 46% He. 


We see that a straight line is obtained. From our 
measurements we obtain for the thermal diffusion 
factor, « = 0.149. 

In Fig. 6 and 7 are shown graphically as an 
example of our measurements the change of the 
velocity as a function of time at 90.1 and 77.3°K 
respectively. 

As the thermal diffusion for H,-He mixtures 
is very small it is difficult to decide from the 
observed decrease whether this is exponential as 
CLARK JoNES and Furry’s calculations indicate. 


Fig. 6. Velocity as a function of time at 90.1°K in H,—He 
mixtures. 
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Fig. 7. Velocity as a function of time at 77.3°K in H,—He 
mixtures. 
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In order to check this theory we reconsidered 
again the measurements carried out by one of us 
together with Forrez and Marrens [1] on 
H.-N, mixtures, for which the thermal diffu- 
sion is much larger than for H,-He mixtures. 

From these curves we calculate the variation 
of concentration as a function of time (see Fig. 9). 


33.3% No; 66.7% He 
P=21.0 cm Hg 


100 200 


300 min 400 


Time 


Fig. 8. Change of the velocity W as a function of time in 
H.-N, mixtures. 
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Fig. 9. Change of concentration as a function of time in 
H.-N, mixtures. 


Corresponding with CLarK Jones and Furry’s 
theory [4] we must have 


(Ay)t = (Ay)1= 00 (1 — e-4) 


with b=CD,,, C being a constant depending 
on the length and the cross-section of the 
connecting tube and also on the temperatures 
T, and T,. 

As the ordinary diffusion coefficient D,, is 
inversely proportional to the pressure p, b 
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must also change in this manner. In Fig. 10 we 
have 
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as a function of t. We observe that straight lines 
are obtained. From these straight lines we 
calculate the values of 6. Thus we obtain: 


b56.8 cm: bs1.0 et 10 Lelie 


which shows b to be proportional to p”. 

We take the opportunity to express our sincere 
thanks to the Belgian Ministry of Education for 
his financial aid during these measurements. 


(Received November 27th, 1954.) 
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GASEOUS OSCILLATIONS 
WITHIN A PERTURBED SPHERICAL RESONATOR 


by I. D. CAmMPBELL* 


Imperial College, South Kensington 


Summary 

A perturbation method is used to investigate theoretically the modes of oscillation of a gas 
within an almost spherical resonator. It is shown that the pressure distribution may be repre- 
sented by the pure mode together with small contributions from other “forced” modes. Under 
some circumstances the resonant frequencies are unperturbed, to the order of approximation 
adopted. 


Sommaire 

On utilise une méthode de perturbations pour étudier théoriquement les modes d’oscillation 
d’un gaz dans un résonateur presque sphérique. On montre que la distribution de la pression peut 
étre représentée par le mode pur auquel s’ajoutent des petites contributions provenant d’autres 
modes «forcés». Dans certaines conditions, les fréquences de résonance ne sont pas perturbées, 
au degré d’approximation considéré pres. 


Zusammenfassung 

Mit Hilfe einer Stérungsmethode wurden die Schwingungsformen eines Gases in einem an- 
nihernd kugelférmigen Resonator theoretisch untersucht. Es wird gezeigt, daf} die Druck- 
verteilung durch die Verteilung einer Eigenschwingung mit kleinen Beitragen anderer ,,er- 
zwungener* Schwingungen dargestellt werden kann. Unter gewissen Umstinden sind die Resonanz- 


frequenzen in der GréBenordnung der hier verwendeten Naherung ungestort. 


1. Introduction 


The modes of oscillation of a loss-free gas 
within a closed spherical resonator have been 
discussed by Ray.ercu [1] and, more recently, 
by Ferris [2]. They show that the pressure 


* Scientific Officer of the Commonwealth of Australia, 
Department of Supply, now at Underwater Research La- 
boratory, Auckland, New Zealand. 


distribution p for a gas oscillating with a pul- 
satance w in a sphere of radius a may be expressed 


p =Jjikr) Pi(u) ett (1) 


where r is the radial co-ordinate, k the wave 
number, and yw is cos) where @ is the polar angle. 
Also j, is the spherical Bessel function and P the 
Legendre polynomial, each of order 1. The 


as 
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resonant frequency of any mode is obtained by 
equating to zero the normal particle velocity at 
the surface of the sphere leading to the equation 


ii(ka) =0. (2) 


The solution of (2) gives a set of values of k, 
and hence w, for each order of mode, /, chosen. 

An experimental investigation of these modes 
led the author to consider the effects of slight 
imperfections in the form of a resonator, nomi- 
nally spherical. Since a search of the literature 
yielded no information the following analysis 
was undertaken. 


2. General analysis 


Consider a deviation from the spherical form 
that is symmetrical about the polar axis. Let the 
surface of the resonator be represented by 


ry =a[l + eg(6)] (3) 


where ro is the radial distance of the surface 
from the origin, a the radius of the sphere, now 
perturbed, ¢ a small dimensionless quantity and 
g(0) a function determining the perturbation. 

If the usual orthogonal spherical coordinate 
system is set up at that point of the surface 
considered, it is found that the normal to the 
surface has the direction cosines 


(1, —edg/dd, 0) (4) 


an expression correct to the first order in e. 
Since the vector particle velocity 0 of an acoustic 
wave is related to the pressure p by 


t= —— Vp (5) 


where @ is the density, and the other symbols are 
as above, the normal particle velocity v, at the 
surface of the resonator is 


i i op € Op dg 
OO or r 06 dé 


(6) 


using (4). The surface of the resonator will be 
assumed rigid, and hence vy, is zero. 

The modes of oscillation of the gas will 
resemble those within a true sphere. The devi- 
ation of any actual mode from the unperturbed 
mode can be represented as the sum of contri- 
butions from other “‘forced”” modes. Hence, for 
a mode symmetrical about the polar axis, let 


P =Jn(s) Pa(u) + € Yam Jnl) P(t) (7) 


where for convenience s has been written for 
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kr. Both here and in the following the time 
factor ej#t has been neglected. Also, without loss 
of generality ~, may be put zero. 

It is convenient to consider a perturbation 
g(9) represented by a Legendre polynomial, say 
P;(), since any arbitrary symmetrical deviation 
may be compounded from a sum of such fune- 
tions. Then the vanishing of vy, at the wall, is 


given by (6) as 
n(s) P(t) — ~ jnls) Pi(u) P¥(u) + 


+ € Lam Jm(s) Pm(u) = 0 


(8) 


where a dash denotes differentiation with respect 
to s, and a star differentiation with respect to 0. 
The variable s is to be evaluated at r= 7p, i.e. at 
the wall. Using (3) to relate ro to a, one obtains 


. Jn(So) P(t) + € 80 Jn(So) Pr() Pi) — (9) 
~ <- In(s0) Pii(u) PR(1u) + & Xam jn( $0) Prn(u) = 0 


where so is ka. The solution of this equation for 
so and a» gives respectively the resonant fre- 
quencies and the coefficients of the forced modes. 
One notes that if ¢ is zero then (9) reduces to (3) 
as expected. 


3. The radial mode 


The solution of (9) is somewhat involved in the 
general case and it is easiest to illustrate the 
solution with particular examples. 

First consider the case of a perturbed radial 
mode, i.e. n is zero. Then (9) becomes (10) 


Jo(so) + € $0 J0(S0) Pr(u) + € ¥%mJm(So) Pm(u) = 0. 


The solution is obtained by equating to zero the 
coefficients of like orders of Legendre poly- 


nominal. For the zero order one has 


(11) 


since ao is zero by supposition. Comparing (11) 
with (2) it is apparent that small deviations from 
a spherical form of resonator leave the resonant 
frequencies of the radial modes unperturbed. 

Similarly, equating to zero the coefficients of 
P)(u) one obtains 


Jo(so) = 90 


So Jo(So) + 1J;(so) = 0 (12) 


where so is the solution of (11). Using the proper- 
ties of the spherical Bessel functions, (12) may 
be transformed to give 


x1 = sojo($o)[j1(s0)- (13) 


It is apparent that all other «,, are zero. So the 
pressure distribution is, ef. (7), 
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F Ss 3 
P= ols) + 050 FH juls) Paw) 
Hence the effect, on a radial mode, of a deviation 
of the shape of the resonator represented by a 
Legendre polynominal of order | is the intro- 
duction of a “‘forced”’ mode of order 1, at the 
unperturbed resonant frequency. 


(14) 


4, The first order perturbation 


Now consider the effect of a perturbation 
represented by the first order Legendre poly- 
nominal on a mode of order n. The resonator is 
now somewhat “egg-shaped”. Equation (9) 
involves the products P, P; and P; Py. It may 
be shown that 


1 
pee ee 0) Pa P,) (15) 
and that 
Tn bh) 
PEPE=— (yi) Pov —Paak (1) 


After substituting (15) and (16) into (9), and 


equating to zero the coefficient of P, one obtains 


Jn(So) = 9 (17) 


showing that the resonant frequenciers are un- 
perturbed. By equating to zero the coefficients of 
P,+41 and P,_; the equations 


(n+1) — jn(so) E 

2n+1) 7441(S0) 

% Me n Jni(S0) | i +1) 
Aad home 1). J),-1(50) 50 


result. Hence the pressure distribution is given 
by 


P =Jnls) Pr(u) + € &n+4Jn41(s) Pn +1(“) + (20) 


— =| , (18) 


Xn+1 = ( 


| (19) 


+ € %n—1Jn—1(8) Pn a(t) 


where x» 1 and xp»_ are as in (18) and (19). If n 
is chosen zero (20) checks with (14). Thus a 
deviation of the shape of the resonator represented 
by a first order Legendre polynomial introduces 
two adjacent forced modes and leaves the reso- 
nant frequencies unperturbed. 


5. The second order perturbation 


Finally consider the effect of a deviation from 
the spherical form represented by the second 
order Legendre polynomial. The shape of the 
resonator now somewhat resembles an hourglass. 

Equation (9) in this case involves the products 
P, Po and Pi P3. In the appendix it is shown 
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that such products may be expanded as series of 
Legendre polynomial by general theorems that 


state that 
n+lI1 
PrPi= ), BaPg 


q=|n—l 


(21) 


and 
n+l 


PPT % vaPa 


fe || 


(22) 


where q takes even or odd values as (n+ l) is 
even or odd. So, in the present case 


Pr Ps = Bn+2 Pre Rie Bn Py i Bn—2 Pr-2 
and 


Pr PS = Yn+2 Pre + Yn Pp + Vn—2 Pn_-2 : (24) 


(23) 


Here it has been assumed that n > 2. The case 
of n= 0 has already been considered, and when 
n=1 the following analysis holds provided f_, 
and y_; are made to vanish. 

Eq. (9) now involves the three polynomials 
Png, Py and Py, 2. The coefficient of P, when 
equated to zero gives 


Jn($0) + € Bn(S0) Jn(So) — & yn (25) 


jul) _ 

$0 
showing that the resonant frequencies are 
perturbed in this case. If the unperturbed values 
of a root of (25) is denoted by s; then the incre- 
mental change 6s; due to the perturbation is 
found to be 


6s, = eae lb, S] + ce 1 

| $1 

Equating to zero the coefficients of Py: one 
obtains 


on +2 = {[1— 


n(n+1)|-) 
a 


a) 


1 
5 sifns2 += 7nsa| x 
HT 


S] 
_jn(si) (27) 
Jn+2(S1) | 


In an analogous way an expression for «,~2 is 
obtained which is (27) with on+2, Bni2, yn+2> 
and jn+2 rewritten as xn-2, Bn—2, Yn-2, and jn—2. 

Thus the pressure distribution of a mode of 
order n in a spherical resonator, the shape of 
which is perturbed by a deviation represented by 
a second order Legendre polynomial is 


P = jn(s) P,(u) +é on +2 Jn+2(s) Pr+2(u) + 
+ €%n—2Jn—2(s) Pn—2(1)- (28) 


6. Discussion 


The three examples considered above illus- 
trate several general conclusions. If the deviation 
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of the shape of the resonator from the spherical 
form is represented by an odd Legendre poly- 
nomial then the resonant frequencies are un- 
changed, to the first approximation. This follows 
from (21) and (22) where no term P,, can occur 
if lis odd. P, can only occur if ] is even and then 
only if n >1/2. So if the surface of the resonator 
is too “corrugated,” (i.e. if 1 is too large), the 
resonant frequency of the mode considered is 
unaffected. 

Secondly, the effect of a deviation represented 
by a Legendre polynomial modifies the “pure” 
mode by introducing other modes forced to 
vibrate at, or near, the resonant frequency of the 
pure mode. If / is odd the number of these forced 
modes is (1+ 1) if n =/ and (n+1) ifn<l. If 
lis even then the number of forced modes is L if 
n=l, nit U25n< land (n--1) atin V2; 

The only case not considered is the trivial one 
of a simple expansion, represented by a zero 
order Legendre polynomial. The change in fre- 
queuncy is easily found from the roots of (2). 

If the deviation g(6), of equation (3), is of 
some general form, symmetrical about the polar 
axis, it is known that it may be represented by a 
series of Legendre polynomials. This series may 
be substituted in (6) and the analysis developed 
in a similar manner to the above. In general the 
resonant frequencies are perturbed and forced 
modes introduced. 


Appendix 


An explicit expression for the coefficients fy 
of (21) was given by Apams [3] and by NEUMANN 
in the same year. The latter’s derivation is 
reproduced by Hopson [4]. In the notation of 
the binomial coefficient 


Roe (get) te 
1 (n+l+ q+) (29) 
n—I+q\ /n+Il—q\ /l—n+q 
(»—1+1)(nit-a) (tne) 
J 2 2 2 


n-+l+g 
fae 
2 


Since Pj, Pf is a polynomial of order (n + 1) in 
4 it is apparent that this product may also be 
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represented as a sum of Legendre polynomials. 
For one has 


dPy, die 
xD ke 78: n l 
PePP==) S244. GO) 
Also 
dP. n(n+1) | 
— 2 UA SS a2 
(1—p?) da (2n +1) (Pn+i1—Pn-1) (31). 
and 
dP 
Ag 7 CID Pi-1+ @l—5) Prat 
3 Py 
+ Q1—9)Prst-t {9 62) 
0 


as | is even or_odd, cf. [4]. The combination of 
(30), (31) and (32) with (21) and (29) leads to the 
expression (22) giving the yg explicitly. 

Alternatively, the approach detailed by 
Hosson [4], (for P, P;) may be adapted to the 
present case. One finds that the product P% Pf 
may be written as 


Ph Pt = ans1(Pn+1+ 02 Pn+i-2+ 


+ w2 04 Pagi-a +++: + 01204..0205 Patios +°*'} 


(33) 


down to a last term in P\,_7;. Explicitly 


vor=—mi(2) G)/PRLF) eo 


an 
Mn+1-28+2—2 
= 35 
se Mn+1-2s ( ) 
where 
M,= (36) 


(r—n—l)(r+n+1+2)(r—n+1+]1)(r+n—I+4 1) 
(2r+1)(r(r+1)—n(n4+1)—l(14+1)} 


; 
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